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ABSTRACT 


The  computer  program  MULTI  accepts  basic  field  inventory  data 
describing  individual  sample  trees.  It  calculates  gross  board-foot  and 
cubic-foot  volume,  grades  or  classifies  for  a  number  of  specified  pri- 
mary products,  adjusts  gross  volume  for  visual  defect,  and  calculates 
standard  error.  Output  tables  indicate  adjusted  gross  volume  per 
acre,  by  grade  and  size  class,  for  each  product  independently.  An 
option  allows  allocation  of  volume  between  several  products  of  a 
multiproduct  combination  following  a  specified  order  of  preference. 

MULTI  was  developed  on  a  CDC  6400  computer  system.  It  is 
written  in  FORTRAN  Extended  Language,  and  can  use  tape  and/or 
disk  storage.  Although  the  program  was  developed  for  products 
and  grading  systems  common  to  ponderosa  pine  (Pinus  ponderosa 
Laws.),  it  is  adaptable  to  other  species,  products,  and  grading  systems. 

Key  Words:  Multiproduct  timber  inventory,  primary  product  evaluation 
(timber),  forest  surveys,  product  volume  estimates,  prod- 
uct volume  adjustments,  timber  utilization. 
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Determining  Timber  Conversion  Alternatives 
Through  Computer  Analysis 


Jack  D.  Heidt,  Donald  A.  Jameson, 
Roland  L.   Barger,  and  Bernard  J.  Erickson 


Timber  Product  Evaluation 

To  achieve  the  best  physical  and  economic 
utilization  of  his  timber  resource,  the  manager 
of  a  diversified  timber  operation  must  be  able 
to  evaluate  the  potential  of  standing  timber 
for  a  variety  of  primary  products.  He  must 
also  be  able  to  estimate  the  effects  of  one 
utilization  program  upon  resource  availability 
or  suitability  for  any  other  program  or  product 
mix.  In  the  long  run,  the  processor  can  try 
to  establish  the  utilization  program  that  makes 
best  use  of  the  available  raw  material.  In  the 
short  run,  however,  he  must  generally  choose 
between  existing  products  and  markets. 

A  method  of  evaluating  standing  timber 
for  a  variety  of  primary  products  has  been 
described  by  Barger  and  Ffolliott  (1970).  Their 
system  may  be  used  with  any  valid  sampling 
scheme,  and  involves  observing  and  recording 
for  each  sample  tree  the  basic  physical  char- 
acteristics that  are  quality  criteria  for  most 
primary  timber  products.  Occurrence  and  sever- 
ity of  defects  form  the  basis  for  estimating  the 
potential  of  the  timber  stand  to  produce  a 
variety  of  products. 

Equally  as  important  as  the  inventory  in- 
formation itself  is  the  ability  of  the  manager 
to  evaluate  these  inventory  data  for  the  entire 
array  of  potential  timber  products.  Large  vol- 
umes of  data  must  be  analyzed  to  determine 
suitability  for  individual  products  and  com- 
binations of  products  in  the  array  of  possible 
conversion  alternatives  (fig.  1).  To  make  prac- 
tical use  of  inventory  data,  the  timber  manager 
or  processor  must  turn  to  automatic  data  pro- 
cessing and  computer  analysis.  The  computer 
program  MULTI  has  been  written  as  a  guide- 
line to  show  one  possible  way  to  implement 
such  a  system. 


MULTI— A  Program  to  Analyze  Product 
Potential 

MULTI  was  developed  on  a  CDC  6400  com- 
puter system.  It  is  written  in  FORTRAN  Ex- 
tended Language,  and  can  use  tape  and/or 
disk  storage. 

The  program  accepts  basic  field  inventory 
data  describing  individual  trees.  It  then  cal- 
culates gross  volume,  grades  or  classifies  for 
a  number  of  primary  products,  adjusts  volume 
for  visual  defect,  and  calculates  standard  error. 

Gross  volume  in  cubic  feet  or  board  feet, 
or  both,  may  be  calculated  for  each  half-log 
or  8-foot  stem  section,  to  minimum  merchant- 
able saw  log  diameter.  Board-foot  volume 
may  also  be  determined  for  16-foot  saw  logs. 
In  addition,  cubic  volume  may  be  calculated 
for  the  stem  section  from  minimum  saw  log 
top  to  a  4.0-inch  top  diameter  inside  bark 
(d.i.b.). 

Depending  upon  the  needs  of  the  user, 
the  program  can  be  used  to  grade  or  classify 
trees  or  sections  of  trees  for  all  possible  prod- 
ucts, or  for  the  individual  products  of  greatest 
value  in  a  specified  order  of  preference. 


Program  Adaptability 

A  key  objective  in  writing  the  computer 
program  MULTI  was  to  develop  a  system  that 
could  easily  be  changed  to  fit  individual  needs. 
Although  the  program  was  developed  for  prod- 
ucts and  grading  systems  common  toponderosa 
pine  (Pinus  ponderosa  Laws.),  it  is  adaptable 
to  other  species,  products,  and  grading  systems. 

Input  parameters  needed  to  calculate  tree 
and  log  volumes  and  estimate  product  volumes 


-  1  - 


Pulpwood 
or 

fiberwood 


Commercial 
poles 

classes  (ft.) 
15 
20 

25-85 


Veneer  logs 

grades 
I 

2 
3 


Saw  logs 

grades 
I 

2 
3 
4 


Stud  logs 

grades 
I 

2 
3 


Combinations 
of  products 

F^guAe  1.--TheJie  oah  many  pKAjnaxy  pKodact  and  pfLoducX-combA^nation  con\j(ifu>lon  a£J:^n.nailvu 
{^oH.  any  poJiticixloA  6tand  o{j  tJjnbdfL.    llZiU,th.aJ:td  oAz  the  pnJjnaJiy  timben.  ptodacJu,  con- 
iid^A^d  in  pn.OQh.am  UULTJ . 


per  acre  are  contained  on  program  data  cards 
and  control  cards.  Individual  products  are 
evaluated  and  graded  through  a  series  of  product 
evaluation  subroutines.  Different  species  or 
products  can  be  accommodated  by  substituting 
appropriate  input  parameters  and  product  evalu- 
ation subroutines.  Changes  in  grading  systems 
require  only  that  the  grading  section  of  the 
affected  subroutine  be  changed. 

For  The  Non-Programer — 
A  Program  Description 

The  following  information  will  describe  for 
the  non-programer  what  data  are  used  and 
what  is  accomplished  by  the  program. 


Field  Data  Input 

Field  data  recorded  for  each  sample  tree 
include: 


1.  Descriptive  and  size  data,  including 

a.  Species. 

b.  Diameter  breast  high  (d.b.h.). 

c.  Total  tree  height  in  merchantable  half  logs. 

d.  Pole  height  for  potential  commercial  poles. 

2.  Visual  defect  characteristics,  including 

a.  Stem  form   defects— sweep,  crook,  fork, 
lean. 

b.  Stem  scar  defects— basal,  lightning. 

c.  Knot  or  limbing  characteristics  in  the 
first  32  feet. 

How  these  data  are  measured  and  recorded  is 
described  in  detail  in  appendix  A. 

Product,  Grade,  and  Priority  Specification 

The  program  can  accommodate  the  full 
range  of  primary  products  and  product  grades 
appropriate  for  a  particular  timber  type  or 
species.  The  program  illustrated  includes  and 
is  limited  to  five  primary  products  from  ponder- 
osa  pine:    saw  logs,  stud  logs,  veneer  logs. 
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commercial  poles,  and  pulpwood.  Within  these 
products,  the  program  will  accommodate  four 
saw-log  grades,  three  stud-log  grades,  three 
veneer-log  grades  (plus  an  "unacceptable"  cate- 
gory in  veneer  and  stud  logs),  and  15  pole 
height  classes. 

Product  priorities  can  be  specified  by  means 
of  a  program  control  card.  If  each  tree  or 
section  is  to  be  graded  and  classified  for  as 
many  products  as  it  is  suitable  for,  no  priority 
selection  is  used.  The  output  will  reflect  the 
total  potential  of  the  stand  for  each  product 
as  if  the  entire  stand  were  being  utilized  for 
that  product.  More  often,  however,  a  descending 
order  of  product  priority  will  be  specified.  A 
typical  order  might  be  (a)  commercial  poles; 
(b)  saw  logs  grade  1,  2,  3;  (c)  veneer  logs  grade 
1,  2,  3;  (d)  stud  logs  grade  1,  2,  3;  (e)  saw  log 
grade  5;  (f)  pulpwood.  In  the  program  illus- 
trated, pulpwood  always  assumes  lowest  priority. 
Each  tree  or  section  would  then  be  graded  for 
these  products  in  order  of  preference  stated, 
and  placed  only  in  the  highest  product  class 
for  which  it  qualifies.  If  desired,  a  series 
of  different  product  priorities  can  be  specified, 
to  be  considered  one  after  another. 


Volume  Determination 

Gross  cubic-foot  and  board-foot  volume  of 
each  sample  tree  is  calculated  from  volume 
equations  appropriate  for  the  species.  Board- 
foot  volume  is  determined  to  minimum  mer- 
chantable saw-log  diameter  in  the  tree,  while 
cubic  volume  is  determined  to  a  4-inch  top 
diameter.  Average  volume  distribution  values 
are  then  used  to  allocate  volume  by  log  and 
half-log  stem  sections.  Volume  equations  and 
volume  distribution  values  used  for  southwest- 
ern ponderosa  pine  are  shown  in  appendix  A. 

Volume  suitable  for  log  products  is  normally 
expressed  in  board  feet,  while  volume  of  pulp- 
wood is  expressed  in  cubic  feet.  Either  ex- 
pression can  be  used  with  any  product,  how- 
ever. For  products  normally  considered  by 
number  or  count,  such  as  commercial  poles, 
a  piece  count  is  maintained. 


Product  Grading  and  Volume  Adjustment 

Each  sample  tree  is  graded  for  the  products 
specified  by  using  a  combination  of  recorded 
field  data  and  program  grading  instructions. 
Grading  instructions  are  based  upon  grading 


or  quality  specifications  for  each  product  (ap- 
pendix A),  and  utilize  the  measured  stem  defect 
characteristics  that  form  these  grading  criteria. 

Tree  size  or  estimated  stem  section  diameter 
is  used  to  judge  initial  suitability  for  each 
primary  product.  Trees  or  stem  sections  of 
acceptable  size  are  graded  for  the  products 
specified.  Gross  volume  graded  is  then  ad- 
justed for  defects  that  either  prohibit  use  of 
the  stem  or  stem  section  for  the  product,  or 
reduce  merchantable  volume. 


Commercial  poles.  —All  sample  trees  9.0 
through  20.9  inches  in  diameter,  for  which  a 
pole  height  has  been  field  recorded,  are  con- 
sidered for  commercial  poles.  The  programed 
evaluation  procedure  for  a  potential  commercial 
pole  proceeds  as  follows: 

1.  Verify  acceptable  stem  size  limits— 9.0  through 
20.9  inches  d.b.h. 

2.  Accept  recorded  pole  height  from  input  field 
data. 

3.  Divide  pole  height  by  8  to  find  8-foot  stem 
sections  or  half-logs  involved.  All  fractions 
of  sections  are  rounded  upward  (pole 
height/8  +  0.9)  to  assure  that  all  sections 
wholly  or  partially  contained  in  the  pole 
are  included. 

4.  Scan  defect  data  for  included  stem  sections 
for  inadmissible  or  limiting  defects. 

a.  The  stem  is  dropped  from  further  con- 
sideration as  a  pole  by  the  occurrence  of 
—sweep  (class  1  or  class  2) 
—lightning  scar 

—lean  of  10°  or  more  (class  2  or  3) 
—knots  4  inches   or  larger  in  diameter 
—major  crook  (class  2)  or  fork  in  any 

except  the  first  or  last  stem  section 

included  in  the  pole  height 
—indication  of  rot. 

b.  Reduce  merchantable  pole  height  by  one 
class  (5  feet)  if  recorded  defects  include 
—fire  or  basal  scar 

— fork  or  major  crook  (class  2)  in  the 
first  or  last  stem  section  included  in 
the  pole  height. 

5.  Pole  count  is  accumulated  by  tree  diameter 
class  and  pole  height  class. 

Pole  height  is  normally  recorded  in  the 
field  only  for  stems  that  do  not  have  inad- 
missible defect.  Programed  defect  screening 
is  consequently  a  "double  check"  on  the  entry 
made  in  the  field.  The  pole  height  recorded 
in  the  field  will  usually  be  accepted  as  correct 
pole  potential  for  the  stem. 
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Saw  logs.— All  sample  trees  of  sawtimber 
size  (11.0-inch  d.b.h.  and  larger  in  ponderosa 
pine)  are  evaluated  for  saw  log  potential.  Saw 
log  grades  are  determined  for  lower  logs,  to  a 
maximum  of  two  logs  or  32  feet  of  merchant- 
able height.  Grades  are  arbitrarily  assigned 
to  logs  above  that  height.  Gross  saw  log 
volume  for  each  log  in  the  tree  is  adjusted  for 
the  presence  of  scaling  defects.  The  programed 
evaluation  proceeds  as  follows: 

1.  Verify  acceptable  stem  size— 11.0-inch  d.b.h. 
or  larger,  plus  merchantable  height  of  two 
or  more  half-logs. 

2.  Determine  gross  board-foot  volume  of  the 
tree,  and  volume  by  16-foot  saw  logs  (in- 
cluding a  top  half-log  where  present),  as 
previously  described. 

3.  Estimate  the  grade  of  each  saw  log  in  the 
tree.  Procedures  for  estimating  grade  depend 
upon  the  position  of  the  log  in  the  tree  and 
the  tree  species  or  grading  system  used. 
Ponderosa  pine  saw  logs  are  graded  as  follows: 

a.  First  16-foot  saw  log: 

Scan  recorded  knot  data  by  8-foot  section 
and  face,   starting  with  stem  section  1, 
face  1,  and  proceeding  in  turn  to  faces 
2,  3,  and  4.     For  each  face,  accumulate 
a  "clear  4-foot  panel"  count  of: 
—no  knots  =  count  two  (2) 
—one  knot  =  count  one  (1) 
—two  or  more  knots  =  count  zero  (0) 
Repeat  for  stem  section  2,  accumulating 
the  panel  count  for  the  first  two  stem  sec- 
tions  (totaling  16  feet  of  merchantable 
length).      From  the  accumulated  panel 
count,  assign  a  log  grade  to  the  log  as 
described  in  appendix  A. 

b.  Second  16-foot  saw  log: 

Grade  the  second  full  saw  log  in  a  tree 
by  the  same  process.  Add  the  clear  panel 
count  for  stem  section  3  to  panel  count 
for  stem  section  4,  and  assign  a  grade. 
Panel  count  for  stem  section  4  is  limited 
to  counting  2  for  each  clear  face  (coded  0), 
since  knot  numbers  are  not  recorded. 
Where  the  third  stem  section  is  the  last 
merchantable  saw  log  volume  in  the  tree 
(that  is,  merchantable  tree  height  of  three 
half-logs),  the  panel  count  for  the  one 
section  is  multiplied  by  two  to  determine 
grade. 

c.  Third  log  and  above: 

Assign  lowest  grade  (grade  5)  to  all  saw 
logs  above  32  feet  in  the  tree.  (Alter- 
natively, upper  saw  logs  can  be  carried 
as  ungraded  volume,  or  can  be  assigned 
other  grades  based  on  available  data.) 

4.  Adjust  gross  volume  of  each  saw  log  for 
visual  defect.    Defects  considered  to  reduce 


usable  volume,  and  average  scale  deductions 
applied,  are: 

Scaling  deduction 
Defect  16-foot  8-foot 

log         (top)  log 

(percent) 


a.  Lightning  scar 

(class  1)  in  tree  25 

b.  Lightning  scar 

(class  2)  in  tree  50 

c.  Sweep  (class  2)  in  tree  20 

d.  Fire  scar  (class  2) 

in  tree  13 
(reduce  butt  log 
scale  only) 

e.  Crook  (any  class)  in 

either  section  of  log  25 

f .  Fork  in  either  section 

of  log  25 


25 

50 
20 


50 


50 


5.  Accumulate  both  gross  and  adjusted  saw 
log  volumes,  by  tree  diameter  class  and  log 
grade. 


Veneer  logs.— All  8-foot  stem  sections  with 
an  estimated  diameter  of  10  inches  or  more 
are  potentially  suitable  for  veneer.  Veneer 
log  quality  evaluation  is  limited  to  the  first 
three  stem  sections  (24  feet),  however,  since 
knot  data  needed  for  grading  are  limited  to 
these  sections.  Veneer  log  evaluation  proceeds 
as  follows: 

1.  Verify  acceptable  stem  section  diameter.  In 
the  absence  of  actual  scaling  diameters,  stem 
sections  in  the  first  24  feet  considered  to 
have  scaling  diameters  of  10  inches  or  more 
are: 


Tree  d.b.h.  class 

12 
14 

16  or  over 


Acceptable  stem  sections 

1st  section 

1st,  2nd  sections 

1st,  2nd,  3rd  sections 


2.  Determine  gross  board-foot  volume  of  the 
tree,  and  volume  of  acceptable  8-foot  stem 
sections,  as  previously  described.  (Merchant- 
able stem  volume  above  24  feet  may  also  be 
calculated  and  carried  as  ungraded  potential 
veneer  log  volume,  if  so  desired.) 

3.  Estimate  the  grade  of  each  potential  veneer 
log  in  the  tree.  The  grade  of  each  8-foot 
log  (maximum  of  three  logs)  is  estimated 
by  scanning  recorded  knot  type  and  size 
data  for  the  four  faces.  Using  size  of  largest 
dead  and  largest  live  knot,  assign  a  grade 
as  described  in  appendix  A. 
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4.  Adjust  gross  volume  of  each  veneer  log  for 
visual  defect.  Defects  considered  to  reduce 
usable  volume,  and  average  scale  deductions 
applied,  are: 

Defect  Scaling  deduction 

(percent) 


a.  Indication  of  rot  in  section  100 
(If  by  d.b.h.  core,  apply 
reduction  to  butt  section  only) 

b.  Crook  (class  2)  in  section  100 

c.  Crook  (class  1)  in  section  50 

d.  Fork  in  section  50 

e.  Lightning  scar  (class  2)  in  tree  50 

f.  Lightning  scar  (class  1)  in  tree  25 

g.  Fire  scar  (class  2)  in  tree  25 
(Apply  reduction  to  butt 

section  only) 


Scan  for  occurrence  of  defects  in  the  order 
shown.  If  total  deduction  reaches  100  per- 
cent, omit  further  scanning  and  apply  100 
percent  deduction. 
5.  Accumulate  both  gross  and  adjusted  veneer 
log  volumes  by  tree  diameter  class  and  log 
grade. 


Stud  logs.— Stud  logs  are  typically  processed 
thi'ough  chipping  headrig  mills,  or  mills  similarly 
equipped  for  the  specialty  manufacture  of  studs 
and  dimension  stock.  The  logs  are  often  pro- 
cessed in  8-foot  lengths,  and  are  restricted 
in  maximum  diameter  by  the  nature  of  the 
processing  equipment.  All  8-foot  stem  sections 
with  an  estimated  scaling  diameter  of  17  inches 
or  less  are  potentially  suitable  for  stud  logs. 
Stud  log  quality  evaluation  is  limited  to  the 
first  three  stem  sections  (24  feet),  however, 
because  of  the  limitations  in  grading.  Stud 
log  evaluation  proceeds  as  follows: 

1.  Verify  acceptable  stem  section  diameter.  In 
the  absence  of  actual  scaling  diameters,  stem 
sections  in  the  first  24  feet  considered  to 
have  scaling  diameters  of  17  inches  or 
less  are: 


by  scanning  recorded  knot  count  and  size 
data  for  the  four  faces.  Using  size  of  largest 
dead  and  largest  live  knot,  and  knot  count, 
assign  a  grade  to  the  stud  log  as  described 
in  appendix  A. 
4.  Adjust  gross  volume  of  each  stud  log  for 
visual  defect.  Defects  considered  to  reduce 
usable  volume,  and  average  scale  deductions 
applied,  are: 

Defect  Scaling  deduction 


(percent) 

a.  Indication  of  rot  in  section  100 
(If  by  d.b.h.  core,  apply  reduction 

to  butt  section  only) 

b.  Crook  (class  2)  in  section  100 

c.  Crook  (class  1)  in  section  50 

d.  Fork  in  section  50 

e.  Lightning  scar  (class  2)  in  tree  50 

f.  Lightning  scar  (class  1)  in  tree  25 

g.  Fire  scar  (class  2)  in  tree  25 
(Apply  reduction  to  butt 

section  only) 


Scan  for  occurrence  of  defects  in  the  order 
shown.  If  total  deduction  reaches  100  per- 
cent, omit  further  scanning  and  apply  100 
percent  deduction. 
5.  Accumulate  both  gross  and  adjusted  stud 
log  volumes,  by  tree  diameter  class  and 
log  grade. 


Pulpwood. — The  cubic  volume  of  all  trees, 
to  a  4-inch  top  diameter,  is  potentially  suit- 
able for  pulpwood.  Pulpwood  is  a  low-value 
product,  however,  and  is  usually  obtained  from 
trees  or  stem  sections  not  suitable  for  other 
primary  products.  Pulpwood  evaluation  is  con- 
sequently limited  to  cubic  volume  not  accepted 
for  other  products,  including  top  volume  from 
minimum  saw  log  top  to  a  4-inch  diameter, 
plus  any  tree  diameter  classes  or  stem  sections 
specifically  excluded  from  other  products.  Pulp- 
wood evaluation  includes: 


Tree  d.b.h.  class 

18  or  under 

20 

22 

24  or  over 


Acceptable  stem  sections 

1st,  2nd,  3rd  sections 
2nd,  3rd  sections 
3rd  section 
None 


2.  Determine  gross  board-foot  volume  of  the 
tree,  and  volume  of  acceptable  8-foot  stem 
sections,  as  previously  described. 

3.  Estimate  the  grade  of  each  potential  stud 
log  in  the  tree.  The  grade  of  each  8-foot 
log  (maximum  of  three  logs)  is  estimated 


1.  Specify  tree  diameter  classes  or  stem  sections 
(if  any)  excluded  from  prior  product  con- 
sideration: smaller  tree  diameter  classes 
for  which  knot  data  are  not  recorded,  tops 
to  4-inch  d.i.b. 

2.  Determine  gross  cubic-foot  volume  of  tree 
or  stem  sections  designated  for  pulp.  For 
sample  trees  11.0  inches  d.b.h.  or  over  which 
have  been  evaluated  for  log  products,  only 
top  volume  will  be  required. 

3.  Adjust  gross  volume  of  each  tree  or  section 
for  visual  defect.     Defects  considered  to 
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reduce  usable  volume,  and  average  volume 
deductions  applied,  are: 


Defect 


Volume 
deduction 

(percent) 

a.  Lightning  scar  (any  class)  in  tree  100 
(including  top) 

b.  Fire  scar  (any  class)  in  tree  50 
(Applies  only  as  reduction  in 

butt  section  volume  when  entire 
tree  is  allocated  to  pulp) 

Scan  for  occurrence  of  defects  in  the  order 
shown.  If  total  deduction  reaches  100  percent, 
omit  further  scanning  and  apply  100  percent 
deduction. 

4.  Accumulate  both  gross  and  adjusted  pulp- 
wood  volumes,  by  tree  diameter  class. 


Other  products.— Potential  for  other  primary 
conversion  products  can  also  be  evaluated  as 
other  products  become  of  interest.  An  analysis 
of  suitability  and  quality  can  be  made  for  any 
product  for  which  the  recorded  stem  charac- 
teristics are  grading  criteria.  Potential  for  a 
specific  product  can  also  be  reevaluated  after 
a  grade  change  or  revision  of  grading  rules. 


Form  of  the  Output 

The  program  prints  out  the  product  analyses 
in  a  series  of  tables.  Each  table  presents 
potential  by  tree  diameter  class  for  one  product 
and  one  or  more  grades  within  the  product 
(appendix  B).  Average  gross  or  adjusted 
volume  per  acre  and  standard  error  are  printed 
out  for  each  grade  within  each  diameter  class, 
and  for  the  accumulated  grade  totals.  Errors 
for  totals  are  generally  sufficient  for  most 
purposes. 

As  an  integral  part  of  the  program,  a  per- 
manent unit  record  is  also  developed  for  each 
sample  tree.  The  unit  record  consists  of  all 
descriptive  variables  recorded  in  the  field,  and 
variables  computed  in  the  course  of  product 
evaluations  (appendix  A).  Thereafter,  output 
information  can  be  obtained  from  the  unit 
records  for  sample  trees,  classified  or  grouped 
by  any  recorded  variables  desired.  Examples 
might  be  (1)  volume  of  grade  3  and  better 
saw  logs  in  blackjack  trees  (as  opposed  to 
old-growth  trees);  (2)  volume  of  acceptable 
veneer  logs  in  trees  with  grade  5  saw  logs 
in  lower  sections;  (3)  volume  of  acceptable 
veneer  logs  in  trees  acceptable  as  commercial 
poles. 


For  the  Programer — 
Program  Implementation 


Calculation  of  Errors  of  Estimate 

The  adjusted  gross  volume  of  each  stem 
section  or  tree  graded  is  converted  to  a  volume- 
per-acre  estimate  by  means  of  the  expansion 
or  conversion  factors  appropriate  for  the  sample. 
For  all  products,  adjusted  gross  volume  or  piece 
count  per  acre  is  then  accumulated  by  product, 
grade,  and  tree  diameter  class,  for  each  random 
sampling  unit. 

Volumes  per  acre  for  each  random  sampling 
unit  provide  the  basis  for  calculating  precision 
of  estimate  for  each  product.  The  program 
uses  a  simple  algorithm  to  compute  standard 
error  of  the  mean  for  each  product  and  grade, 
by  tree  diameter  class.  The  algorithm  is  most 
useful  for  simple  random  sampling  and  sub- 
sampling  with  units  of  equal  size.  It  can  also 
be  used  for  proportional  stratified  sampling. 
Subsampling  systems  with  substantially  unequal 
sample  units  may  require  weighting  procedures. 

Calculated  error  terms  are  ultimately  printed 
out  with  the  volume  estimates  to  which  they 
correspond. 


Program  MULTI  operates  as  an  executive 
routine  that  provides  for  initial  data  card  read-in, 
tree  diameter  class  computation,  initial  tests  of 
sample  tree  diameter  against  product  size  limit 
criteria,  and  preparation  of  a  tape  or  disk  file 
of  tree  records  for  further  processing  (fig.  2). 
Once  the  initial  data  have  been  tape  or  disk 
filed,  MULTI  is  designed  to  make  multiple 
passes  through  the  data,  executing  options  as 
specified  on  program  control  cards.  MULTI 
calls  on  various  subroutines  to  do  the  neces- 
sary grading,  volume  computations,  and  eiTor 
calculations.  An  option  is  provided  for  writing 
a  unit  record  tape  file  for  each  tree  processed. 

The  program  utilizes  eight  subroutines.  Five 
are  the  product  evaluation  routines  POLE,  SAW, 
VENR,  STUD,  and  PULP,  with  the  names 
corresponding  to  the  product  being  evaluated. 
The  remaining  three  are  VOLUM,  which  cal- 
culates board-foot  and  cubic-foot  volumes  of 
sample  trees,  CVOL,  which  converts  sample 
volumes  to  volume  per  acre,  and  STDER,  which 
calculates  confidence  limits  around  the  means. 
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A  complete  printout  of  program  MULTI,  in- 
cluding all  subroutines,  is  included  as 
appendix  B. 


Program  Data  Cards 

Program  data  cards  type  A,  B,  and  C  are 
read  in  at  the  beginning  of  the  program.  These 
cards  provide  the  board-foot  and  cubic-foot 
volume  distribution  percentages  needed  to  cal- 
culate volume  by  stem  section  (see  appendix  A). 

Card  type  A  includes  a  set  of  three  cards 
containing  cubic-foot  volume  distribution  per- 
centages for  half-log  (8-foot)  stem  sections, 
plus  top  to  4-inch  d.i.b.  The  data  are  used  to 
fill  the  array  KCUV8(65).  Card  type  B  includes 
three  cards  containing  board-foot  volume  dis- 
tribution percentages,  also  by  half-log  (8-foot) 
stem  section.  The  data  are  used  to  fill  the 
array  KBFV8(55).  Card  type  C  includes  two 
cards  containing  board-foot  volume  distribution 
percentages  for  16-foot  saw  logs.  Data  are  used 
to  fill  the  array  KBFV16(29).  In  each  case, 
volume  distribution  data  are  included  for  tree 
heights  and  log  or  half-log  positions  through 
10  half-logs  or  five  saw  logs. 

Complete  program  card  data  for  southwest- 
ern ponderosa  pine  are  shown  in  appendix  A. 
The  cards  are  coded  as  follows  (format  2613), 
with  all  values  right-adjusted  in  the  data  fields: 


Card 

Card 

Card 

Percent 

Array 

type 

number 

columns 

value 

element 

A 

A-1 

1-3 

65 

L(l,l) 

4-6 

35 

L(1,T) 

7-9 

50 

L(2,l) 

etc. 

etc. 

etc. 

73-75 

10 

L(6,5) 

76-78 

8 

L(6,6) 

A-2 

1-3 

3 

L(6,T) 

4-6 

26 

L(7,l) 

7-9 

18 

L(7,2) 

etc. 

etc. 

etc. 

73-75 

8 

L(9,7) 

76-78 

6 

L(9,8) 

A-3 

1-3 

4 

L(9,9) 

4-6 

2 

L(9,T) 

B 


B-1 


B-2 


B-3 


C  C-1 


C-2 


7-9 

19 

L(10,l) 

etc. 

etc. 

etc. 
. 

34-36 

4 

L(10,10) 

37-39 

2 

L(10,T) 

1-3 

100 

L(l,l) 

4-6 

66 

L(2  1) 

7-9 

34 

L(2,2) 

10-12 

52 

L(3,l) 

etc. 

etc. 

etc. 

73-75 

13 

L(7,4) 

76-78 

10 

L(7,5) 

1-3 

7 

L(7,6) 

4-6 

6 

L(7,7) 

7-9 

25 

L(8,l) 

etc. 

etc. 

etc. 

73-75 

9 

L(10,6) 

76-78 

7 

L(10,7) 

1-3 

6 

L(10,8) 

4-6 

4 

L(10,9) 

7-9 

4 

L(10,10) 

1-3 

100 

L(l,l) 

4-6 

75 

L(1.5,l) 

7-9 

25 

L(1.5a.5) 

10-12 

64 

L(2,l) 

etc. 

etc. 

etc. 

73-75 

35 

L(5,l) 

76,78 

27 

L(5,2) 

1-3 

20 

L(5,3) 

4-6 

13 

L(5,4) 

7-9 

5 

L(5,5) 

The  volume  distribution  data,  illustrated 
here  for  ponderosa  pine,  are  input  by  means 
of  cards  to  allow  the  user  to  select  the  set  of 
distribution  percentages  that  he  wishes  to  use. 
If  values  for  other  species  are  desired,  they  can 
be  substituted  in  the  same  manner.  If  addi- 
tional values  are  required  which  exceed  the 
dimensions  of  the  arrays  mentioned,  array  size 
should  be  increased  and  the  source  code  modi- 
fied or  expanded  in  initial  READ  statements 
and  in  subroutine  VOLUM.  The  data  cards 
type  A,  B,  and  C  are  initial  inpvit  only,  and 
are  not  repeated  for  successive  passes. 
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Program  Control  Cards 

Program  control  and  option  parameters  are 
read  in  on  four  control  cards  (card  types  1 
through  4)  at  the  beginning  of  each  problem 
and  before  each  pass  through  the  data. 
Program  control  cards  are  coded  as  follows, 
with  all  values  right-adjusted  in  the  data  fields: 


Variable  Card 
name  columns 

PROB  1-80 


IN 


1-5 


NSAM  6-10 


R 


11-15 


NPROD  16-20 


IP(I)  21-75 
EXPFC  76-80 


Card  Type  1 


Format  Description 

8A10  Alphameric  code  to 
describe  problem 
and  user  options. 

Card  Type  2 

15       IN  =  0,  read  tree 
data  from  cards. 
IN  =   1,  read  tree 
data  from  tape. 

15  Number  of  sample 
units.  If  "cluster" 
sampling  is  used, 
NSAM  is  number  of 
clusters. 
F5.3  If  R  is  punched  to 
any  non-zero  value, 
rot  occurrence  ob- 
served and  recorded 
in  a  subsample  of 
trees  will  be  ran- 
domly extended  to 
the  entire  sample. 
If  rot  occurrence 
data  are  available 
for  all  sample  trees, 
R  should  be  punch- 
ed zero. 

15  Number  of  products 
to  be  evaluated  in 
each  run.  Maximum 
size  of  NPROD  is 
11,  1  pole  grade,  4 
saw  log  grades,  3 
veneer  grades,  and 
3  stud  grades. 
1115  Sets  the  product  pri- 
ority order  for  pro- 
duct evaluation. 
F5.0  EXPFC,  when  equal 
to  some  value  other 
than  zero,  is  a  con- 
version factor  to 
convert  fixed-size 
plot  sample  values 
to  per-acre  basis.  In 
plotless    or  point 


Variable  Card 

name    columns  Format  Description 

sampling,  computed 
conversion  factors 
are     used,  and 
EXPFC  =  0. 
Format  (215,  F5.3,  1215,  F5.0) 


JON  1-5 


KPOLE  6-10 


KSAW  11-15 


KVENR  16-20 


KSTUD  21-25 


KPULP  26-30 


KOUT  31-35 


IBUG  36-40 


Card  Type  3 

15  JON  =  1,  specifies 
product  evaluation 
on  priority  basis, 
each  stem  section 
evaluated  only  for 
highest  priority 
product  for  which 
it  is  suited,  follow- 
ing specified  order 
of  preference  per 
IP(I),  card  type  2. 
JON  =  2,  specifies 
primary  product  e- 
valuation  in  which 
each  stem  section  is 
evaluated  for  all 
products  for  which 
it  is  suited. 
JON  =  3,  specifies 
both  evaluations  to 
be  made. 

15  Trees  to  be  e- 
valuated  for  poles, 
KPOLE  =  1;  other- 
wise KPOLE   =  0. 

15  Trees  to  be  e- 
valuated  for  saw 
logs,  KSAW  =  1; 
otherwiseKSAW=0. 

15  Trees  to  be  e- 
valuated  for  veneer, 
KVENR  =  1;  other- 
wise KVENR  =  0. 

15  Trees  to  be  evalu- 
ated for  stud  logs, 
KSTUD  =  1;  other- 
KSTUD  =  0. 

15  Trees  to  be  evalu- 
ated for  pulpwood, 
KPULP  =  1;  other- 
wise KPULP  =  0. 

15  KOUT  =  1,  write  on 
output  unit  record 
tape;  KOUT  =  0, 
skip  this  option.  A 
unit  record  can  be 
written  only  if 
JON  =  2  or  3. 

15  IBUG  =  1,  print  out 
individual  tree  sta- 
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tistics;  IBUG  =  0, 
skip  this  option. 
Use  only  for  testing 
or  debugging  pro- 
gram as  a  large  a- 
mount  of  output  is 
generated. 

KVOL     41-45  15       KVOL  =  1,  compute 

and  output  gross 
volume  tables; 
KVOL  =  2,  compute 
and  output  adjusted 
gross  volume  tables, 
volumes  adjustedfor 
visual  scaling  de- 
fects. 

BAFV     48-55        F8.4       BAF  value  of  the 

prism  or  angle  gage 
used  in  plotless 
sampling. 

Format  (915,  2X,  F8.4) 

Sample  input  for  a  production  run  might  be: 
3-1-1-1-1-1-1-0-1-25.0000 

Card  Type  4 

POINT(I)1-80  16F5.0  Specifies  the  num- 
ber of  sampling 
points  in  each  sam- 
ple unit  (NSAM).  If 
more  than  16  sample 
units  are  used,  con- 
tinue on  subsequent 
cards.  If  cluster 
sampling  is  not 
used,  POINT  (I) 
should  equal  1  for 
all  entries.  Decimal 
point  need  not  be 
punched. 

Product  Priority  Instructions 

If  products  are  to  be  evaluated  on  a  priority 
basis,  the  desired  order  of  priority  must  be 
specified  by  entries  in  the  IP(I)  fields,  columns 
21-75,  of  control  card  type  2.  Entries  refer 
to  the  specific  subroutine  being  called: 

1  Subroutine  POLE 

2  Subroutine  SAW 

3  Subroutine  VENR 

4  Subroutine  STUD 

If  the  desired  priority  is  poles,  followed  by 
grades  1-3  saw  logs,  grades  1-3  veneer  logs, 
grades  1-3  stud  logs,  and  grade  5  saw  logs, 
the  IP(I)  values  would  be  coded  as  follows: 

Variable  Column  Number  Product  considered 

NPROD  19-20  11 

IP(1)  25  1  Poles 

IP(2)  30  2  Saw  logs,  grade  1 


ir(o) 

oD 

9 

9 

k_J Cl  VV    1             ^  b 

o 

V  tJllwCl  ,  gld.U.t;  L 

xr  } 

o 

IP(7) 

55 

3 

Veneer,  grade  3 

IP(8) 

60 

4 

Studs,  grade  1 

IP(9) 

65 

4 

Studs,  grade  2 

IP(IO) 

70 

4 

Studs,  grade  3 

IP(ll) 

75 

2 

Saw  logs,  grade  5 

Computed  Control  Variables 

A  number  of  control  variables  are  com- 
puted during  the  execution  of  the  program. 
These  variables  are  as  follows: 

KMAX  Equals  the  maximum  diameter 
class  in  each  data  set. 

KMIN  Equals  the  minimum  diameter 
class  in  each  data  set. 

IND(I)        Is  an  indicator  which  is  set  to 

1=1,  10  one  to  indicate  that  a  particular 
log  or  tree  has  been  utilized 
for  a  priority  product.  IflTRHT 
is  greater  than  10,  ITRHT=10. 

KOUNT(I)  Equals  the  number  of  trees  in 
each  diameter  class. 

NBORE  Number  of  trees  sampled  for 
rot  determination. 

NROT  Number  of  sampled  trees  con- 
taining rot. 

PROROT  Proportion  of  trees  with  rot  as 
computed  by  NROT/NBORE. 
The  program  uses  a  random 
number  generator  to  randomly 
extend  the  observed  frequency 
of  occurrence  of  rot  to  the  entire 
sample.  If  rot  occurrence  data 
are  available  for  all  sample  trees, 
the  program  provides  an  option 
(R=0  )  to  bypass  this  procedure 
and  use  rot  data  directly.  In 
either  case,  occurrence  of  rot 
should  be  coded  2  for  proper 
program  evaluation. 


Tree  Data  Cards 

Inventory  data  are  recorded  on  punch  cards, 
one  card  for  each  sample  tree.  The  data  in- 
cluded on  the  cards,  and  the  coding  and  format 
required  for  the  program,  are  described  and 
illustrated  in  appendix  A. 

As  many  tree  data  cards  as  desired,  within 
the  limits  of  the  disk  or  tape  file,  may  be  read 
in.  Data  may  be  entered  in  any  order  without 
regard  to  diameter  class  or  sample  number. 
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One  blank  card  is  used  to  signal  the  end  of 
tree  card  input;  two  additional  blank  cards 
(read  as  type  1  and  2  control  cards)  will  termi- 
nate the  program  after  the  desired  computation 
and  output  is  completed. 

If  additional  runs  or  problems  are  to  be 
executed  on  the  filed  data,  repeat  program 
control  cards  types  1,  2,  and  3.  Control  card 
type  4  is  not  required  for  additional  runs. 
As  many  sets  of  type  1,  2,  and  3  control  cards 
may  be  input  as  desired,  with  the  last  set 
being  followed  by  two  blank  cards. 


Subroutine  POLE  (PLHT2,  POINT,  TOTAL, 
ISAM,  IMAX,  PLHTl) 

The  arguments  PLHT2,  POINT,  and  TOTAL 
are  dimensioned  variables,  the  sizes  set  by  the 
argument  ISAM,  which  is  equal  to  or  greater 
than  NSAM,  the  number  of  sample  units  used. 
The  argument  PLHTl  accumulates  POLE  counts 
by  diameter  class.  The  argument  PLHT2  accumu- 
lates POLE  counts  by  sample  number  and  pole 
height  class.  In  other  product  subroutines, 
each  stem  section  is  evaluated  independently 
for  the  product  being  considered,  and  if 
JON  =  1  accepted  sections  are  withdrawn  from 
consideration  for  subsequent  products.  In  sub- 
routine POLE,  however,  the  use  of  a  tree  as  a 
pole  eliminates  any  consideration  for  other 
products  and,  conversely,  use  of  any  8-foot 
section  for  another  product  eliminates  the  entire 
tree  from  being  considered  as  a  pole. 


Subroutine  SAW  ( VCNT,  POINT,  TOTAL,  ISAM, 
IMAX) 

The  arguments  VCNT,  POINT,  and  TOTAL 
are  dimensioned  variables,  the  sizes  set  by 
ISAM  and  IMAX.  ISAM  must  be  equal  to  or 
greater  than  the  total  number  of  samples 
(NSAM),  and  IMAX  must  be  equal  to  or  greater 
than  the  largest  diameter  class  in  the  data. 
The  array  VCNT  is  an  accumulator  for  saw  log 
volumes  by  diameter  class  and  grade,  POINT 
is  a  vector  specifying  the  number  of  points 
(or  subplots)  in  each  sample,  and  TOTAL  is 
an  array  for  accumulation  of  volumes.  This 
subroutine  evaluates  each  sample  tree  between 
a  minimum  diameter  class  (KMIN)  and  maxi- 
mum class  (KMAX)  that  has  at  least  one  saw 
log  (two  8-foot  sections).  Descriptive  tree 
quality  data  are  screened  by  8-foot  stem  section, 
and  a  grade  from  the  designated  grading  system 
is  assigned  each  saw  log.  The  grading  system 
and  method  of  evaluating  grading  criteria  in 
ponderosa  pine  are  described  in  appendix  A. 


If  priorities  are  established  (JON  =  1),  selection 
of  either  or  both  half-logs  for  a  higher  priority 
will  eliminate  the  full  log  from  saw  log  con- 
sideration. Log  volumes  can  be  computed 
and  printed  out  in  terms  of  either  gross  volume 
(KVOL  =  1)  or  volume  adjusted  for  visual 
scaling  defect  (KVOL  =  2),  as  previously  dis- 
cussed. Information  is  printed  out  by  diameter 
class  and  total,  for  one  log  grade  at  a  time 
if  JON  =  1  or  for  all  grades  at  once  if  JON  =  2. 

Subroutine  VENR  (VCNT,  POINT,  TOTAL, 
ISAM,  IMAX) 

Input  arguments  for  this  subroutine  are  the 
same  as  those  described  for  SAW.  This  sub- 
routine evaluates  each  sample  tree  meeting 
initial  diameter  class  specifications;  it  screens 
by  8-foot  stem  section  and  assigns  a  grade 
from  the  designated  grading  system  (appendix 
A).  Since  grades  are  based  on  recorded  knot 
data,  grading  is  limited  to  the  first  three  8-foot 
stem  sections  in  each  tree.  Log  volume  can 
be  adjusted  for  visual  scaling  defects  (KVOL  =2) 
if  desired.  If  priorities  are  established 
(JON  =  1),  printouts  will  include  only  the 
three  acceptable  grades.  A  "grade  4"  (material 
unacceptable  because  of  knot  size)  is  accumu- 
lated, however,  and  under  option  J0N=2  will 
be  printed  out  and  listed  in  the  tree 
unit  record. 


Subroutine  STUD  (VCNT,  POINT,  TOTAL, 
ISAM,  IMAX) 

The  input  arguments  and  operations  of  this 
subroutine  are  the  same  as  those  described  for 
VENR.  The  same  options  and  restrictions  apply 
to  grade  evaluation  (appendix  A)  and  print- 
outs for  the  product. 


Subroutine  PULP  (PULPV,  POINT,  TOTAL, 
ISAM,  IMAX) 

The  argument  PULPV  is  a  dimensioned 
variable  of  size  ISAM.  The  other  arguments 
are  the  same  as  in  other  subroutines.  PULP 
accumulates  cubic-foot  volumes  for  all  sample 
tree  sections  not  accepted  for  other  products, 
plus  top  volume  from  minimum  saw  log  top 
to  a  4-inch  top  diameter  for  each  tree  evaluated. 
All  trees  excluded  because  of  d.b.h.  less  than 
9  inches  are  also  included  in  pulpwood  totals. 
Pulpwood  gross  volumes  can  be  adjusted  for 
inadmissible  visual  defect  (fire  and  lightning 
scar)  as  previously  discussed,  under  the 
KVOL  =  2  option. 
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Subroutine  VOLUM  (KCUV8,  KBFV8,  KBFV16, 
BAFV) 

This  subroutine  calculates  total  merchant- 
able board-foot  and  cubic-foot  volume  in  each 
sample  tree.  Board-foot  volume  is  calculated 
to  minimum  saw-log  diameter  and  cubic-foot 
volume  is  calculated  to  a  4-inch  top  diameter. 
Volumes  are  allocated  by  log  and  half-log  stem 
sections,  by  means  of  average  volume  distribu- 
tion values.  Tree  volume  equations  and  volume 
distribution  percentages  used  for  ponderosa 
pine  are  described  in  appendix  A.  The  argu- 
ments KCUV8,  KBFV8,  and  KBFV16  are  di- 
mensioned variables  (volume  distribution  per- 
centage values)  read  in  by  the  program  data 
cards. 

If  plotless  or  point  sampling  has  been  used, 
the  subroutine  also  calculates  expansion  factors 
to  convert  sample  values  to  per-acre  values. 
An  expansion  factor  is  calculated  for  each  sample 
tree,  from  the  equation 

F.F  =  BAFV 

0.005454D' 

I where  BAFV  is  input  on  control  card  type 
3  and       is  tree  d.b.h.  squared. 


Subroutine  CVOL  (VBL,  EXPFC,  EF) 

This  is  a  function  subroutine  used  to  con- 
vert product  count  or  volume  per  tree  or  stem 
section  to  count  or  volume  per  acre,  when 
plotless  sampling  has  been  used.  Argument 
VBL  is  a  dummy  argument  used  to  pass  the 
sample  value  to  be  converted  to  the  function 
CVOL.  Argument  EF  is  the  point  sample 
per-acre  expansion  factor  calculated  by  sub- 
routine VOLUM. 

EXPFC  (also  called  PRISM  in  some  sub- 
routines) is  a  factor  to  convert  fixed-size  plot 
sample  values  to  per-acre  values.  EXPFC  =  0 
when  plotless  sampling  is  used. 


Subroutine  STDER  (NSAM,  Y,  S,  AMTRX, 
POINT) 

This  subroutine  uses  a  simple  algorithm  to 
compute  standard  error  of  the  mean  for  each 
product  and  grade,  as  indicated  earlier.  Cal- 
culated error  terms  are  printed  out  with  the 
volume  estimates  for  the  product.  Arguments 
for  the  subroutine  are: 

Input 

NSAM     Number  of  sample  units 
AMTRX  Data  vector  of  size  NSAM 
POINT    Vector  of  size  NSAM  giving  number 
of  points  or  subplots  per  sample. 
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Return 

Y         Mean  of  vector  AMTRX 

S         Standard  error  of  vector  AMTRX. 


Unit  Record  Output 

A  permanent  unit  record,  containing  descrip- 
tive variables  and  a  number  of  computed  vari- 
ables, is  developed  for  each  sample  tree.  The 
variables  included  in  the  record  are  described 
in  appendix  A.  Unit  record  output  is  gen- 
erated when  program  control  card  type  3  speci- 
fies JON  =  2  or  3  and  KOUT  =  1.  The  record 
is  written  on  external  device  logical  file  4 
(tape  or  disk  storage)  by  a  WRITE(4)  LIST 
statement,  and  consists  of  strings  of  binary 
word  values  in  the  form  in  which  they  appear 
in  storage.  Each  tree  record  consists  of  105 
binary  words.  An  end-of-file  is  written  after 
the  last  tree  processed;  a  message  is  then 
printed  that  gives  the  number  of  records  in 
the  output  file,  plus  a  printout  of  the  last 
tree  record  processed  for  check  purposes. 

It  is  assumed  that  the  unit  record  output 
file  would  normally  be  used  as  input  on  the 
CDC  6400  system  for  further  processing  by  the 
IOCS  input-output  routines.  If  this  is  not 
the  case,  it  would  be  advisable  to  write  the 
unit  record  output  with  a  WRITE(4,FMT)  LIST 
statement  in  the  BCD  mode.  BCD  tapes  written 
with  this  statement  can  normally  be  read  on 
a  different  computer,  often  with  the  normal 
FORTRAN  input-output  routines. 


Program  Operation  and  Limitations 

Program  MULTI  is  written  in  the  FORTRAN 
Extended  Language  Version  1.0,  an  extension  of 
the  USASI  FORTRAN  language,  for  the  Control 
Data  6400  computer  system.  Operation  of  the 
program  requires  a  central  processor,  a  card 
reader,  output  line  printer  (assigned  logical 
file  3)  and  two  external  tape  or  disk  storage 
devices  (assigned  logical  files  1  and  4). 

The  program  contains  error  messages  which 
provide  information  to  help  the  user  correct 
the  input  data.  There  are  no  recovery  points  in 
the  program.  When  an  abort  or  execution 
failure  occurs,  all  volume  accumulations  to 
time  of  failure  are  lost  and  the  program  must 
be  restarted  from  the  beginning. 

Maximum  diameter  class  and  sample  size 
are  set  in  the  data  statement,  DATA  IMAX, 
ISAM/40.4,.  Dimensioned  variables  SAWV, 
VENRV,  STUDV,  PULPV,  VCNT,  PLHT2, 
TOTAL,  and  POINT  contain  as  their  first  sub- 


script  the  maximum  sample  size  allowed. 
Because  all  of  these  variables  except  VCNT 
are  intermeshed  with  an  EQUIVALENCE  state- 
ment, their  first  dimensions  must  be  the  same 
and  should  equal  ISAM.  The  last  dimension 
of  SAWV,  VENRV,  and  VCNT  must  equal  IMAX. 
Any  change  in  ISAM  and/or  IMAX  requires 
a  program  coding  change  in  all  the  variables 
mentioned.  Other  dimensioned  variables  do 
not  require  changes  due  to  changes  in  sample 
number  or  maximum  diameter  class. 

A  test  version  of  progi^am  MULTI,  with 
maximum  sample  size  =  130  and  maximum 
diameter  class  =  34,  required  a  field  length 
of  105,7008  0"^"  core  storage  of  36, o  K  to 
load  the  program.  A  CDC  6400  computer  sys- 
tem with  140,000  8  or  50, o  K  available  can 
possibly  accommodate  maximum  variable  com- 
binations to  sample  size  =  150  and  diameter 
class  =  40.  Larger  combinations  may  require 
overlay  subroutines. 

Loading  procedure 

The  order  for  loading  the  program  is: 

1.  Computer  center  control  cards. 

2.  MULTI-Main  and  all  subroutines. 

3.  Program  data  cards,  types  A,    B,    and  C 
(a  set  of  eight  cards) 

4.  Program  control  cards,  types  1,  2,  3,  and  4. 

5.  Individual  tree  data  cards. 

6.  Blank  card  to  signal  end  of  tree  data. 

7.  Two  additional  blank  cards  if  run  is  to  termi- 
nate after  first  pass,  or 

Repeat  program  control  cards,  types  1,  2,  and 


3  as  many  times  as  desired,  one  set  for  each 
pass. 

8.  Two  blank  cards  to  terminate  job,  final  pass. 
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APPENDIX  A 


Field  Inventory  Data 


Field  inventory  methods  are  described  in 
Barger  and  Ffolliott  (1970).  Field  inventory 
data  described  here  are  those  taken  in  inven- 
torying ponderosa  pine.  For  other  species  or 
primary  products,  some  modification  of  the  data 
recorded  may  be  desirable.  For  the  evaluation 
program  MULTI,  one  data  card  is  punched 
for  each  sample  tree,  following  the  coding 
and  format  described: 


Variable  Card 
name  columns 

—  1-2 

ISTRT  3-4 

IDMY  5-11 

DBH  12-16 

IPLHT  18-19 


Format  Description 


ISWP 

IPCR 

ICRK 
IPFK 
IFRK 
IFS 

ILS 


KNOT 

(1,1,1) 
KNOT 
(1,1,2) 


21 

22 

23 
24 
25 
26 

27 

29 
30 


2X 
12 
17 

F5.2 
12 

II 
II 

II 

II 
II 
II 

II 

II 

II 


Blank  or  identifica- 
tion data. 
Sample  number, 
range  01  to  99. 
Dummy  variable, 
can  be  used  for  ID 
data. 

Tree  d.b.h.,  punch 
decimal  point, 
XX. XX. 

Pole  height,  record- 
ed as  maximum  5- 
f  oot  height  class  ob- 
tainable (25',  30', 
35',  etc.). 

Sweep  entered  as  1 
(minor)  or  2 
(major). 

Number  of  half-log 
in  which  crook 
occurs. 

Crook  entered  as  1 
(minor)  or2  (major). 
Number  of  half-log 
in  which  fork  occurs. 
Coded  as  1  if  fork 
exists. 

Fire  scar  entered  as 
1  (minor)  or  2 
(major). 

Lightning  scar  en- 
tered as  1  (minor) 
or  2  (major). 
Number  of  knots  in 
face  1,  section  1. 
Size  of  largest  live 
knot  to  nearest  inch, 
face  1,  section  1. 


Variable  Card 
name  columns 

KNOT  31 
(1,1,3) 

KNOT  32-34 
(1,2,N) 

KNOT  35-37 
(1,3,N) 

KNOT  38-40 
(1,4,N) 

KNOT  43-54 
(2,M,N) 


KNOT  57-68 
(3,M,N) 


I24P 

I24P(1)  71 

I24P(2)  72 

I24P(3)  73 

I24P(4)  74 

ITRHT  76-77 


LEAN  79 


Format  Description 

II      Size  of  largest  dead 
knot  to  nearest  inch, 
face  1,  section  1. 
311       Repeat  for  face  2, 
section  1,  N=l,2,3. 
311       Repeat  for  face  3, 
section  1,  N=l,2,3. 
311      Repeat  for  face  4, 
section  1,  N=l,2,3. 

1211  Same  knot  data  re- 
corded by  8-foot  sec- 
tion and  face  for  sec- 
ond section  (or  to 
minimum  merchant- 
able top  diameter). 

1211  Same  knot  data  re- 
corded by  8-foot  sec- 
tion and  face  for 
third  section  (or  to 
minimum  merchant- 
able top  diameter). 

For  the  fourth  8-foot 
section  (24'  -  32'), 
only  presence  or  ab- 
sence of  clear  8-foot 
faces  is  entered  as: 

II  Face  1:  0  if  clear, 
1  if  knots  occur. 

II  Face  2:  0  if  clear, 
1  if  knots  occur. 

II  Face  3:  0  if  clear, 
1  if  knots  occur. 

11  Face  4:  0  if  clear, 
1  if  knots  occur. 

12  Tree  height,  record- 
ed as  total  half-logs 
or  8-foot  sections 
to  the  minimum 
saw  log  top  diam- 
eter. If  tree  height 
is  recorded  in  to- 
tal feet  or  in 
16-foot  sections,  a 
conversion  to  half- 
log  height  is  re- 
quired. 

II  Lean,  recorded  by 
5  degree  classes 
(nearest  class)  as: 
0  -  Less  than  3  de- 
grees 
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Variable    Card  Format 
name  columns 


Description 


1  -  5  degrees 

2  -  10  degrees 

3  -  15    degrees  or 

more 

IROT         80  II        Rot,  as  observed  in 

subsample  by  means 
of  increment  cores 
or  other  indicators, 
recorded  as: 

0  -  not  sampled 

1  -  sampled,  no  rot 

2  -  sampled,  rot 
Card  format  (2X12,  17,  F5.2,  1x12,  1x711,  1x1211, 
2x1211,  2x1211,  2x411,  1x12,  1x211) 


Sample  Tree  Volume  Determination 

Cubic-foot  Volume  Determination 

Determine  tree  height  in  half-logs  (8-foot 
sections)  to  minimum  merchantable  saw  log 
diameter.    If  only  total  tree  height  has  been 


recorded,  estimate  half-log  height  from  total 
height  (Van  Deusen  1967). 

Calculate  total  merchantable  cubic-foot  vol- 
ume of  sample  tree  to  4.0-inch  top,  by  equation. 
For  tree  height  in  16-foot  logs  to  nearest  half- 
log  (H)  and  d.b.h.  outside  bark  (D),  equations 
for  southwestern  ponderosa  pine  are  (Myers 
1963): 


(1) 


(2) 


Tree  d.b.h.  under  19.0  — 
D^H  800  or  less: 

V  =  0.046000  D^H  +  6.800000 
D  ^  H  over  800  : 


V  =  0.044204  D'H 
Tree  d.b.h.  over  19.0  — 
D^H  1,000  or  less  : 

V  =  0.050666  D 
D^H  over  1,000  : 

V  =  0.045736  D'H  ^ 


8.266000 


H  +  5.866800 


10.857212 


Nineteen  inches  is  the  assumed  breaking  point 
between  blackjack  and  old-growth  trees. 

Determine  volume  distiibution  among  stem 
sections  of  the  tree.  Apply  percentages  to  total 
tree  volume  to  estimate  volume  of  sections. 
Average  cubic  volume  distribution  in  south- 
western ponderosa  pine  of  specified  half-log 
height  is: 


Tree 

Half- 

log 

position  in 

tree 

Top  to 

height 

4 

0- inches 

ha  If -  logs 

1 

2 

3 

4 

5 

6  7 

8 

9 

10 

d.i.b. 

Number 

Percent  of 

total 

cubic 

-foot 

tree 

volume 

I 

65 

35 

2 

50 

31 

19 

3 

43 

29 

17 

11 

4 

38 

27 

17 

11 

7 

5 

32 

24 

17 

13 

9 

5 

6 

28 

21 

16 

14 

10 

8 

3 

7 

26 

18 

15 

13 

10 

9  6 

3 

8 

24 

16 

14 

12 

10 

9  7 

5 

3 

9 

22 

15 

13 

12 

10 

8  8 

6 

4 

2 

10 

19 

15 

12 

11 

10 

8  8 

6 

5 

4 

2 

Board-foot  volume  determination 

Determine  tree  height  in  half-logs  (8-foot 
stem  sections)  to  minimum  merchantable  saw- 
log  diameter.  If  only  total  tree  height  has  been 
recorded,  estimate  half-log  height  from  total 
height. 

Calculate  total  merchantable  board-foot  vol- 
ume of  sample  tree  to  minimum  merchantable 
top  diameter  (variable),  by  equation.  For  tree 
height  in  16-foot  logs  to  nearest  half-log  (H) 
and  d.b.h.  outside  bark  (D),  Scribner  rule  equa- 


tions for  southwestern  ponderosa  pine  are 
(Myers  1963): 

(1)  Tree  d.b.h.  under  19.0— 


(2) 


D  '  H  800  or  less: 


H  +  8.165600 


V  =  0.224793  D 
D  '  H  over  800: 

V  -  0.300081  D  '  H  -  52.090112 
Tree  d.b.h.  19.0  or  over  — 

D  '  H  1,130  or  less: 

V  =  0.275784  D  '  H  -  5.091250 
D  ^  H  over  1,130: 

V  =  0.326427  D  '  H  -  62.962331 
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Nineteen  inches  is  the  assumed  breaking  point 
between  blackjack  and  old-growth  trees. 

Determine  volume  distribution  among  stem 
sections  of  the  tree.  Average  board-foot 
volume  distribution  in  southwestern  ponderosa 
pine  of  specified  log  height  is  indicated  below. 
For  saw  logs,  normally  expressed  in  terms  of 


16-foot  logs  and  half-logs,  refer  to  the  fii'st 
tabulation.  For  veneer  and  stud  logs  or  other 
log  products  utilized  in  8-foot  lengths,  refer 
to  the  second  tabulation.  Apply  percentages 
to  total  tree  volume  to  estimate  volume  of 
individual  logs  or  sections. 


Tree 

height  in 

16-foot 

saw  log  position  in  tree 

Logs 

Half-logs 

1 

1.5  2 

2.5        3      3.5        4  4.5 

5 

Number 

Percent  of 

total  board— foot  tree  volume  - 

1 

2 

100 

1.5 

3 

75 

25 

2 

4 

64 

36 

2.5 

5 

58 

32 

10 

3 

6 

49 

33 

18 

3.5 

7 

41 

32 

19  8 

4 

8 

44 

31 

17        --  8 

4.5 

9 

39 

29 

19        --       10  3 

5 

10 

35 

27 

20        --  13 

5 

Tree 

8-foot 

half- 

log  position  in  tree 

height, 

half- logs 

1 

2 

3 

4 

5 

6 

7      8  9 

10 

Numbe  r 

Percent 

of 

total 

board 

-foot 

tree  volume 

1 

100 

2 

66 

34 

3 

52 

32 

16 

4 

43 

30 

18 

9 

5 

36 

26 

18 

12 

8 

6 

31 

23 

18 

14 

9 

5 

7 

28 

20 

16 

13 

10 

7 

6 

8 

25 

18 

15 

13 

10 

9 

6  4 

9 

23 

16 

14 

12 

10 

9 

7      5  4 

10 

21 

15 

13 

11 

10 

9 

7      6  4 

4 

Product  Specifications  and  Grades 
Ponderosa  Pine  Commercial  Poles 

All  trees  9.0  through  20.9  inches  in  diam- 
eter and  of  acceptable  pole  form  are  con- 
sidered potential  poles.  Grading  specifications 
for  commercial  poles  were  adapted  from  those 
of  the  Ameiican  Standards  .Association  (1963). 

Inadmissible  defects  in  commercial  poles 
include 

—sweep  (deviation  greater  than  1/3  d.b.h.) 
—major  crook  (deviation  greater  than  1/2 

pole  diameter  at  crook) 
—knots  larger  than  4  inches  in  diameter 

(dead  or  green) 


—knot  whorls  or  clusters  aggregating  more 
than  8  inches  of  knot  diameter  within 
1  linear  foot 

—fork 

—heart  rot 

—lightning  scar 

—fire  scar 

—compression  wood 
in  stems  leaning  10' 

For  trees  meeting  minimum  merchantable 
specifications,  stem  length  to  the  first  limiting 
defect,  such  as  length  to  first  inadmissible  knot, 
is  field  recorded.  Such  defects  as  fork,  crook, 
and  fire  scar,  if  located  near  the  butt  or  top  of 
the  stem,  may  not  eliminate  the  pole  but  re- 
quire a  reduction  of  acceptable  pole  length. 


(considered  present 
or  more) 


-17- 


Programed  grading  procedures  scan  defect 
data  for  stem  sections  included  in  pole  height, 
and  either  verify  or  change  recorded  pole  height 
as  necessary. 


Ponderosa  Pine  Saw  Logs 

All  trees  11.0  inches  d.b.h.  and  larger  are 
considered  sawtimber  trees,  and  logs  8.0  inches 


and  larger  in  scaling  diameter  are  evaluated 
as  saw  logs. 

All  grading  specifications  apply  to  16-foot 
log  lengths.  The  same  specifications  apply  to 
shorter  logs,  in  proportion  to  their  length. 

All  logs  meeting  the  minimum  meixhant- 
ability  standards  are  graded  by  the  improved 
grading  system  for  ponderosa  pine  logs  (Gaines 
1962). 

Abbreviated  grading  specifications  are  as 
follows: 


Defects  Permitted 


Grade  Primary  (log  knots) 

1  One  log  knot  not  over  %  inch 

in  diameter. 


Secondary  (scar,  etc.) 

Confined  to  three  4-foot  panels 
or  less. 


2  Confined  to  four  4-foot  panels  Secondary  plus  primary  confined 
or  less.            "  to  six  4-foot  panels. 

3  Six  4-foot  panels  free  of  all  grading  defects. 

4  (Logs  of  the  type  described  by  grade  4  of  this  system  do  not 
generally  occur  in  southwestern  ponderosa  pine;  consequently, 
the  grade  is  omitted  from  consideration.) 

5  All  other  logs  with  net  scale  of  one-third  or  more  of  gross  scale. 


Programed  grading  instructions  estimate 
the  number  of  clear  4-foot  panels  in  each  log, 
based  on  recorded  knot  data.  Accumulated  panel 
count  is  used  to  assign  a  log  grade  to  the  log, 
as  follows: 

Clear  panels,  Log  grade 
1/4  circumference  x  4  feet 

15-16  1 
12-14  2 
5-11  3 
Less  than  5  5 


Ponderosa  Pine  Veneer  Logs 

Logs  10.0  inches  and  larger  in  scaling  diam- 
eter    are  considered  potential  veneer  logs. 

All  grading  specifications  apply  to  8-foot 
log  lengths. 

Logs  that  meet  the  minimum  merchant- 
ability requirements  for  veneer  logs  are  graded 
according  to  the  following  arbitrary  grading 
rules  for  ponderosa  pine  veneer  logs,  wiitten 
for  use  in  multiproduct  inventory  analysis. 

GRADE  I    Veneer  logs  from  which  a  high 
proportion  of  grade  A  and  B 


veneer  can  probably  be  re- 
covered. 

Grade  Specifications 

Dead  knots  -  none  allowed. 
Green  knots  -  allowed  up  to 
2  inches  in  diameter. 
GRADE  II  Veneer  logs  from  which  a  high 
proportion  of  grade  C  veneer 
can  probably  be  recovered. 
Grade  Specifications 

Dead  knots  -  allowed  up  to 
2  inches  in  diameter. 
Green  knots  -  allowed  up  to 

2  inches  in  diameter. 
GRADE  III  Veneer  logs  from  which  a  high 

proportion  of  grade  D  veneer 
can  probably  be  recovered. 
Grade  Specifications 

Dead  knots  -  allowed  up  to 

3  inches  in  diameter. 
Green  knots  -  allowed  with- 
out size  limit. 

UNACCEPTABLE  Logs  identified  as 
"grade  4"  for  computer 
programing  and  output 
purposes. 

Dead  knots  -  4  inches 
or  more  in  diameter. 
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Ponderosa  Pine  Stud  Logs 

Logs  6.0  through  16.9  inches  in  scaling 
diameter  are  considered  potential  stud  logs. 

All  grading  specifications  apply  to  8-foot 
log  lengths. 

Logs  that  meet  basic  size  requirements  for 
stud  logs  are  graded  according  to  the  following 
arbitrary  grading  rules  for  ponderosa  pine  stud 
logs,  written  for  use  in  multiproduct  inventory 
analysis. 

GRADE  I    Stud  logs  from  which  a  high 
proportion  of  SELECT  and  CON- 
STRUCTION grade  studs  can 
probably  be  recovered. 
Grade  Specifications 

Dead  knots  -  allowed  to  1  inch 

in  diameter. 

Green  knots  -  allowed  to  2 
inches  in  diameter. 
Total  number  of  knots  -  can- 
not exceed  16. 
GRADE  II  Stud  logs  from  which  a  high 
proportion  of  STANDARD  grade 
studs  can  probably  be  recovered. 
Grade  Specifications 
Dead  knots  -  allowed  to  2 
inches  in  diameter. 
Green  knots  -  allowed  to  2 
inches  in  diameter. 
Total  number  of  knots  -  can- 
not exceed  32. 
GRADE  III  Stud  logs  from  which  a  high 
proportion  of  UTILITY  and 
ECONOMY    grade    studs  can 
probably  be  recovered. 
Grade  Specifications 
Dead  knots  -  allowed  to  2 
inches  in  diameter. 
Green  knots  -  allowed  to  3 
inches  in  diameter. 
Total  number  of  knots  -  un- 
limited. 

UNACCEPTABLE  Logs  identified  as 
"grade  4"  for  computer  pro- 
graming and  output  pur- 
poses. 

Dead  knots  -  3  inches  or  more 
in  diameter. 

Green  knots  -  4  inches  or 
more  in  diameter. 

Sample  Tree  Unit  Record 

The  unit  record  for  each  sample  tree  in- 
cludes the  following  descriptive  and  computed 
variables  applicable  to  the  tree,  printed  in  the 
following  order: 

Tree  diameter  class  (K) 

Sample  number  (ISTRT) 


Sample  identification  data  (IDMY) 
Tree  d.b.h.  (DBH) 

Commercial    pole    merchantable  height 

(IPLHT) 
Defect  occurrence: 
Sweep  (Class  1,  2:  ISWP) 
Crook     (Location     and     Class     1,  2: 

IPCR,  ICRK) 
Fork  (Location  and  occurrence:  IPFK, 
IFRK) 

Fire  scar  (Class  1,  2:  IFS) 
Lightning  scar  (Class  1,  2:  ILS) 
Knot  data: 

Count,   size,    first  stem  section  KNOT 
(N,L,D) 

Count,  size,  second  stem  section  KNOT 
(N,L,D) 

Count,  size,  third  stem  section  KNOT 
(N,L,D) 

Clear  face  count,  fourth  stem  section 
I24P(N) 

Half-log  merchantable  tree  height  (ITRHT) 
Lean  (Class  1,  2,  3:  LEAN) 
Indication  of  internal  rot  (IROT) 
Pole  selector  indicator  (0  =  No,  1  =  Yes) 
(KP) 

Total  tree  cubic-foot  volume  to  4-inch  top 
(TCFV) 

Total  tree  board-foot  volume  to  saw  log 

top  (TBFV) 
Vol. /Acre  expansion  factor  used  for  tree 

(EF) 

Volume  of  saw  logs  (board  feet),  by  log 
position  and  grade 
BFSW(1,1),  grade  1,  log  position  1 
BFSW(2,2),  grade  1,  log  position  2 
BFSW(1,3),  grade  1,  log  position  3 
BFSW(1,4),  grade  1,  log  position  4 
BFSW(1,5),  grade  1,  log  position  5 
Repeat  for  saw  log  grades  2,  3,  4  (grade 
4=5) 

Volume  of  veneer  logs  (board  feet),  by 
log  position  and  grade 
BFVN(1,1),  grade  1,  log  position  1 
BFVN(1,2),  grade  1,  log  position  2 
BFVN(1,3),  grade  1,  log  position  3 
Repeat   for   veneer   log  grades    2,  3,  4 
(grade  4  =  unacceptable) 

Volume  of  stud  logs  (board  feet),  by  log 
position  and  grade 
BFST(1,1),  grade  1,  log  position  1 
BFST(1,2),  grade  1,  log  position  2 
BFST(1,3),  grade  1,  log  position  3 
Repeat    for    stud    log   grades    2,  3,  4 
(grade  4  =  unacceptable) 

Volume  of  pulpwood  (cubic  feet);  top  piece 
and  other  sections 

CFPL(l),  in  top,  saw  log  to  4-inch  top 
CFPL(2),  all  sections  not  acceptable  for 
other  products 
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APPENDIX  B 


PROGRAM  MULTI  HUL  10 

1 (  1 NPUT,OUTPUI .TAPE5= INPUT, TAPE3  =  0UTPUT, TAPE  1, TAPE4 )  MUL  20 

C  COMPUTER    PROGRAM   TO  CONVERT    MULTIPROOUCT    INVENTORY   DATA    TO  YIELD        MUL  30 

C  BY  PRODUCT  MUL  40 

C                                                             BY  MUL  50 

C                           JACK   D.    HE!OT   AND  DONALD   A.    JAMESON  MUL  60 

C  ROCKY   MOUNTAIN   FOREST    AND  RANGE    EXPERIMENT    STATION  MUL  70 

C                                                                 PROGRAM  CONTROL    CARDS  MUL  80 

C   MUL  90 

C  DATA  CAROISl    TYPE    A  MUL  100 

C   MUL  110 

C  A   SET  OF   THREE  CARDS   TO   INPUT   THE   PERCENT  OF  TOTAL  CUBIC  FOOT  VOL.   MUL  120 

C  BY  8-FT.    STEM   SECTIONS.   SEE   PROG.    DESCRIPTION  FOR  ORDER  AND  FORMAT  MUL  130 

C   MUL  1^0 

C  DATA  CARO(S»    TYPE  B  MUL  150 

C   MUL  160 

C  A   SET   OF    THREE    CARDS   TO    INPUT    THE    PERCENT   OF    TOTAL    BOARD   FOOT    VOL.    MUL  170 

C  pY  8-FT,    STEM   SECTIONS.    SEE   PROG.    DESCRIPTION   FOR  ORDER   AND  FORMAT  MUL  180 

C   MUL  190 

C  DATA  CARO(S)    TYPE    C  MUL  200 

C   MUL  210 

C  A    SET  OF    TWO  CARDS    TO    INPUT    THE    PERCENT   OF    TOTAL    BOARD   FOOT   VOL.        MUL  220 

C  BY  l6-f=T.    STEM  SECT  I  HNS  (  SAWLOGS  I .    SEE   PROG.   DESCRIPTION  FOR  ORDER     MUL  230 

C  AND  FORMAT  MUL  240 

C   MUL  250 

C  CONTROL  CARD  TYPE   1  MUL  260 

C   MUL  270 

C  PROB(I)=PROBLEM    IDENTIFIER  MUL  280 

C   MUL  290 

C  CONTROL  CARD  TYPE   2  MUL  300 

C   MUL  310 

C  IN=0   INDICATES  TREE   DATA  TO  BE    INPUT   FROM  CARDS  MUL  320 

C  1N=1    INDICATES   TREE   DATA   TO   BE    INPUT    FROM  TAPE  MUL  330 

C  NSAM=NUMBER  OF    SAMPLE   UNITS    IF   CLUSTER    SAMPLING    IS   USED,    NSAM    IS        MUL  340 

C  THE    NUMBER   OF    CLUSTERS.  MUL  350 

C  FOR    PROPORTIONAL   HOT    REDUCTION,    R=NON-ZERO  MUL  360 

C  NPRDD=NUMBER    OF    PRODUCTS    TO   BE    EVALUATED  MUL  370 

C  lP(n=PROnuCT    priority  MUL  380 

C  IF   USING   PRIMARY  PRODUCT   SELECTION  CRITERIA    IPIl)    SPECIFIES  ORDER     MUL  390 

C  OF    PRODUCT    PRIORITIES  HUL  400 

C  FOR   PLOTLESS  CRUISING,   EXPFC=0.    FOR  FIXED  PLOT   SIZE  CRUISING,  MUL  410 

C  EXPFC    IS   A  CONVERSION  FACTOR   TO  ADJUST   RESULTS  BASED  ON  PLOT   SIZF     MUL  420 

C   MUL  430 

C  CONTROL    CARD  TYPE    3  MUL  440 

C   MUL  450 

C  JON=l    INDICATES   PRIMARY   PRODUCT    SELECTION   CRITERIA  MUL  460 

C  J0N=2    INDICATES    SAMPLE   WILL    BE    EVALUATED   FOR   ALL    PRODUCTS  MUL  470 

C  J0N=3  GIVES   BOTH   TYPES   OF    EVALUATION  MUL  480 

C  FOR    POLE    STUDY,    KP0LE=1,    TO  OMIT   KPOLE=0  MUL  490 

C  FOR    SAWLOG    STUDY   KSAW=l,    TO  OMIT    KSAW=0  MUL  500 

C  FOR    VENEER   STUDY,    KVENR=1,    TO  OMIT   KVENR=0  MUL  510 

C  FOR    STUD   LOG    STUDY,    KSTU0=1,    TO   OMIT   KSTUD=0  MUL  520 

C  FOR    PULP    STUDY  KPULP=1,T0  OMIT  KPULP=0  MUL  530 

C  FOR    UNIT    RECORD  OUTPUT    T A PE , JON= 2 0R3 , K OUT= 1 , TO   OMIT   KOUT=0  MUL  540 

C  FOR    INDIVIDUAL   TREE    COMPUTATIONS    IBUG=1.T0  OMIT    ieUG=0  mUL  550 

C  FOR   OROSS  VOL   OUTPUT   KVOL=UFOR   ADJUSTED  GROSS  VOL  OUTPUT  KV0L  =  2        muL  560 

C  UNIT    RECORD  OUTPUT    IS    ACCUMULATED   ONLY   FOR    GROSS    VOL,    KV0L=1.  mUL  570 

C  BAFV=PASAL    area  FACTOR   VALUE   OF    THE   PRISM   OR    ANGLE   GAUGE    USED    IN        muL  580 

C  PLOTLESS  CRUISING  MUL  590 

C   MUL  600 

C  CONTROL   CARD   TYPE    4  HUL  610 

C   MUL  620 

C  THIS   CARD  WILL    6F   READ  ONLY    ON    INITIAL    PASS    THROUGH  DATA  muL  630 

C  pniNT(I)    =   NUMBER   OF    PLOTS    PER   CLUSTER  MUL  640 

C  USE  ADDITIONAL  CARDS    IF   NEEDED  HUL  650 

C  IF   CLUSTER   SAMPLING    IS   NOT  USED,    P0INT(I)  =  1   FOR   ALL   ENTRIES  huL  660 

C   MUL  670 

C  COMPUTED  CONTROL  VARIABLES  MUL  680 

C   MUL  690 

C  KMAX=MAXIMUM   DIAMETER  CLASS    IN   DATA  MUL  700 

C  KMIN=MINIMUM   DIAMETER    CLASS    IN   DATA  HUL  710 

C  IN0(l)=PRinR   USE    INDICATOR  HUL  720 

C  KOUNTI  n=NUMBER    OF    TREES    PER    DIAMETER   CLASS  MUL  730 

C  NflORE=    NUMBER  OF    TREES    80RED  MUL  740 

C  NROT    =    NUMBER   OF    BORED    TREES   WITH   RDT  MUL  750 

C   HUL  760 

C                                                             INDIVIDUAL    TREE    DATA  CARDS  HUL  770 

C   HUL  780 

C  THESE    CARDS   WILL    NOT    BE   READ    IF    IN.NE.O  HUL  790 

C  1STRT=SAMPLE    NUMBER  HUL  800 

C  |OMY=OUMMY   VARIABLE, CAN    BE   USED  FOR   TREE    ID    INFORMATION  HUL  810 

C  DBH=D1AMETER    BREAST   HIGH  MUL  820 

C  IPLHT=POLE    HEIGHT  MUL  830 

C  ISWP=SWEEP   CLASS  MUL  840 

C  IPCR  =  CROOK   LOCATION  HUL  850 

C  ICHK=CROOK  CLASS  MUL  860 

C  iPFK=FORK   LOCATION  MUL  870 

C  |FRK=FORK  CLASS  HUL  880 

C  IFS  =  FIRE    SCAR  MUL  890 

C  ILS=LIGHTNING    SCAR  HUL  900 

C  KNOT(K,L,M)=KNOT    NUMBER    AND   SIZE  MUL  910 

C  I?4P(n=CLEAR  FACE   CLASSES  MUL  920 

C  |TRHT=TREE    HEIGHT    IN   HALF    LOGS  MUL  930 

C  LEAN=LEAN  CLASS  MUL  940 

C  IROT=ROT    CLASS  MUL  950 

C  E^JD   TREE    DATA  CARD    INPUT    WITH   ONE    BLANK    CARD  MUL  960 

C   MUL  970 

C  INDICATORS    INTERNAL    IN   THE    PROGRAM  MUL  980 

C  VARIABLE    DIMENSIONS   OF    ALL    DIMENSIONED  VARIABLES   ARE    SHOWN    IN  HUL  990 

C  DATA    STATEMENT,    SEE    SUBROUTINE   DIMENSIONS    FOR    GUIDE  MULIOOO 

C  NO  CHANGES    ARE    REQUIRED    lit    SUBR.    AS  VARIABLE   DIMENSIONS   ARE    USED  MUHOIO 

C  I^flM   EQUAL    TO   OR   GREATER    THAN  NSAM  MUL1020 

C  IMAX    EQUAL    10   OR   GREATER    THAN  KMAX  MULI030 

C  FOR    UNIT   RECORD  OUT PUT , RE QUE  ST   TAPE4    FOR    OUTPUT.    HAS    EOF    MARK    END  MUL1040 

C   MUL1050 

Dl MENS  I  ON   SAHV(4,4,40) ,VENRV{4,4,401, STUDVI 4,4,9) , PL HT  1(15, 201,  HUL 1060 

IPLHT2(4, 15 ) ,T0TAL(4,4 ) ,P0INT(4) ,PULPV( 4) ,KCUV8( 65) ,KBr V8( 55) ,  MUL 10  70 

2KBFV16(29)  MUL1080 

C  MUL 1090 
COMMON  KNOT! 3,4,3) , I24P (4) , IP( 11 ) , INDl 1 0 ) , KOUNT ( 50 1 , Y I  1 5 ) ,5(15).       MUL 1 100 

IPR08( 8) ,PCUV8( II ) ,PBFV8( 10) , P8F V I  6 ( 5 ) , BF SW ( 4, 5 ) , BFVN I  4, 3 ) ,  MUL 1 1  10 

2eFST14,3»,CFPL(2) . NPR OD , K M I N , KM AX , NS AM , K , ISTRT,DBH, IPLHT, I SwP,  MUL 1 120 
3IPCR,ICRK,IPFK,IFRK,IFS,ILS. ITRHT,LEAN, I R OT , N SA W , N VE NR , NS TUO , JON ,      MUL  I  130 

4LAST.TCFV, TBFV, I  BUG , I C ARD, PR  I SM , KP , IDMY,KVOL,EF  MUL  1  140 

C  MUCH  50 
EQUIVALENCE  (SAWVf  1  ,  I  .  I  )  ,  PL  HT2  (  1  ,  1  H  ,  (  SAWV(  1  ,  1,5)  ,  TOTAL!  1.  U  )  ,  MUH  160 
1  ( SAWV( 1  ,  1  ,6)  ,PULPV( 1  I  )  ,    ( SAWV(  1  ,  1  , 7  )  , PO  1  N T (  1 )  )  ,    ( SAWV(  I, 1 ,9) , VENRV  MUL  1  170 

2(1,1.101  I  MUL I  180 

C  HUL1190 

DATAIMflX, ISAM/40, 4/  MUL1200 

KPASS=0  HUL1210 

NHFC=0  MUL1220 

KP=0  MUL1230 

EF=0.0  HUL 1240 

00    100    1=1,4  MUL1250 


DO  100  J=l.3 

BFVN(l,J)=0.0 
100  BFST(I.J)=0.0 

00   102  1=1.4 

DO   102  J=l,5 
102  BFSH(I,J)=0.0 

CFPL( 1 )=0.0 

CFPL(2)=0.0 

 HEAD  VOLUME   PERCENTAGE   TABLES   FOR  CUFT-8 , BRFT-8 ,    AND  BRFT-16. 

 INPUT  CARDS,    TYPES   A, B.C. 

RE AD (5, 106)  ( KCUV8( 1 ) , 1=1 ,65) 
106   FnRMAT(26I3/(2613) ) 

REA0(5 , 106 1  1 KBFV8 ( n , 1=1 ,55 1 

READ(5 , 106)  ( K6FV16I  I  ), 1  =  1,29  1 

REWIND  4 
999   REWIND  I 

KPASS=KPASS  ♦  1 

 READ  CONTROL  CARDS     CARD  TYPES  1,2,3,4 

READ    (5,40)    (PROB( I ) . 1=1.8) 

READ    (5  ,42)    I N . NS AM , R , NPR OD ,  (  IP(  I  )  .  I  =  1  , NPR 00  )  , E X P FC 

IF    (NSAM. GT. ISAM)    WRITE    (3,411  NSAM 

IF    (NSAM. GT. ISAM)    GO  TO  39 

IF    (NSAM.EO.O)   GO  TO  39 

PRISM=EXPFC 

READ (5,43)    JON , KPOLE . KS AH , KVENR , KSTUO, KPULP ,KOUT, I  BUG, KVOL , BAFV 

JONREF= JON 

J0N=AHIN0(JON,l) 

IF(KPASS.GE.2)   GO  TO  2 

READ    (5,44)    ( PO I  NT ( 1) , 1  =  I , NS AH ) 

2  KMAX=0 
KM1N=100 
NROT=0 
NBORE=0 
PROROT=0.0 

DO  3  K=2,50,2 

3  KOUNT(K)=0 

4  IF    ( IN.NE.O)    GO   TO  5 

 READ  TREE   DATA  CARDS.   THE  FORMAT   AND    INPUT  LIST   MAY  BE  CHANGED  TO 

 AGREE  WITH  USERS  COOED  DATA.    ALL  VARIABLES  NOT    INCLUDED   IN  USERS 

 INPUT   LIST   SHOULD  BE   SET   TO  ZERO  FOR  FURTHER  USE   IN  THE  PROGRAM 

RE AO (5, 45  I ISTRT, lOMV.DBH, IPLHT, ISWP, IPCR. ICRK, IPFK, IFRK, IFS. ILS, ( 
1 (KNOT(L,M,N) ,N=1 ,3) ,M=1 ,4},L=1, 3), ( I24P(  I ), 1  =  1, 4 ) , I TRHT , LE AN , I  ROT 
IFl I TRHT.GT.  10)    ITRHT  =  10 

 COMPUTE  DIAMETER  CLASS 

K=OBH/2.0+0.5 
K=K*2 

IF    (K.GT.IMAX)    WRITE    (3,46)  OBH 
IF    (K.GT. IMAX)   GO  TO  4 

IF    ( ISTRT. GT. ISAM)    WRITE    (3,47)  ISTRT 
IF    (  ISTRT. GT. ISAM)   GO   TO  4 


C  RECORD  TREE  DATA  ON  TAPE 


WRITE ( 1 IK,ISTRT, IDMY,DBH,  IPLHT, ISWP,  IPCR,  ICRK, IPFK,  IFRK, IFS, ILS.  I  ( 
1 (KNOT(L,M,N) ,N=1,3) ,M= 1,4) ,L=1, 3 ) , (  I24P( 1  ) , I  =  1 , 4 ) , I TRHT , LE AN t I  ROT 

GO   TO  6 


 IF   DATA   IS  TAPE   FILED,   READ  FROM  TAPE 

5  READ( 1  IK, ISTRT, IDMY.OflH, IPLHT,  ISWP, I  PC R, ICRK,  IPFK, IFRK, IFS, ILS, (!( 
IKNOTCL.M.N) ,N=1,3) ,M=1,4) ,L  =  1,3),(I24P(I),I  =  1,4),I TRHT, LEAN, I  ROT 

6  IF    (D6H)  11,11,7 

7  KOUNT I K 1 =KOUNT ( K ) *1 

IF    (IROT.GT.Ol  NB0RE=NB0RE+1 
IF    (IROT.GT.l)  NRnT=NROT+l 
I F    ( K-KMAX )  9,9,8 

8  KMAX=K 

IF    (K-KMIN)  10,4,4 
10  KMIN=K 

GO  TO  4 
U    REWIND  1 

IF    (R.NE.OI    PRQROT=FLOAT (NROT)/ FLOAT INBORE) 
ICARD=0 

GO  TO   112,26,12),  JONREF 


-BEGIN  PRODUCT    PRIORITY  ROUTINE 


-READ   AND  EVALUATE    DATA   BY  PRODUCTS 


DO  25  KX=KHIN,KMAX,2 

LIM=KOUNT (KX ) 

IF    (LIM.EQ.O)   GO  To  25 

00  24    J=l ,LI M 

DO  13  1=1,10 

13  IND(  n  =0 

READ11)K,ISTRT,I DM Y, OBH,  I PLHT, ISwP,  IPCR,  ICRK.  IPFK,  I FRK , I F S , I L S ,  (  (  ( 
1KN0T( L,M,N) ,N=1.3) ,M=1,4) ,L  =  1, 3) , (  I24P1 I  )  ,  I  =  1  , 4 ) ,  I  TRH T , L E AN , I  ROT 
I F( I TRHT.LE.O)    WRI  TE( 3.47  1  )    I TRHT, IDMY 
IF ( I TRHT.LE.O)   GO  TO  24 

CALL    VOLUM  (KCUV8,KBFVR,KBFV16.flAFV) 

 IF   RANDOM  NUMBER  GE NE R AT OR ( R ANF ( 1 ) )    IS  NOT  USED  TQ  DETERMINE  SUB" 

 SAMPLE   FOR  ROT  REDUCT I  ON , 1 NSERT  A  GO   TO   16   BRANCH  CARD  HERE.  R=0 

GO  TO  16 

IF    (RANF ( 1  )  .LT.PROROT )  14,15 

14  IR0T=2 
GO   TO  16 

15  IRaT=0 

16  IF    ( IBUG.EQ.O)    CO   TO  18 
ICARD=ICAR0*1 

WRI TE( 3,49)    I  CARD, ISTRT, PO I  NT  1  I STRTl ,eF,K,DBH, IPLHT, I TRHT, ISWP, 
1  LEAN, ICRK, IPCR, IFRK, I PFK,  IFS,  ILS,  I ROT, ( 124PI  IM ),  IM  =  1  ,4) 
DO   17  L=l,3 
00    17  M=l,4 

17  WRITE    (3,50)    L , M , ( KNOT ( L , M, N ) , N= 1 , 3 ) 

18  IF    (  IKX.EO.KMAX ) . AND. ( J.EQ.L IM)  )    LAST  =  3 
NSAW=0 

NVENR=0 
NSTUD=0 

DO  23  I=l,NPR00 
ISUB=1 P( I ) 

GO  TO   119,20,21,221,  ISUB 

19  CALL   POLE    ( PLHT2 , POI NT , TOT AL , I S AM ,  I  MAX , PLHT 1 ) 
GO  TO  23 


HUL1260 

MUL1270 

MUL1280 

MUL1290 

MUL 1300 

MUL1310 

MUL 1320 

MUL1330 

MUL1340 

MUL1350 

MUL 1360 

MUL1370 

MUL1380 

MUL  1390 

MUL1400 

MUL1410 

MUL 1420 

MUL1430 

MUL1440 

KUL1450 

MUL1460 

MUL  1470 

MUL  1480 

MUL1490 

MUL1500 

MUL1510 

MUL1520 

MUL 15 30 

MUL1540 

MUL1550 

MUL1560 

MUL1570 

MUl.1580 

MUL1590 

MUL 1600 

MUL1610 

MUL1620 

MUL  1630 

MUL1640 

MUU1650 

HUL  1660 

MUL1670 

MUL1680 

MUL 1690 

MUL1700 

MUL1710 

MUL1720 

MUL 1730 

MUL1740 

MUL1750 

MUL1760 

MUL1770 

MUL 1730 

MUL  1790 

MUL 1800 

MUL  1810 

MUL  1820 

MUL1830 

MUL  1840 

MUL1850 

MUL1860 

MULia70 

MUL1880 

MUL1890 

MUL  1900 

MUL1910 

MUL1920 

MUL1930 

MUL1940 

MUL1950 

MUL I960 

MUL  1970 

MUL1980 

MUL 1990 

MUL2000 

MUL2010 

MUL2020 

MUL2030 

MUL2040 

MUL2050 

MUL2060 

MUL2070 

MUL2080 

MUL2090 

MUL2100 

MUL21 10 

MUL2120 

MUL2130 

MUL2140 

MUL2150 

MUL2160 

MUL2irO 

HUL2180 

MUL2190 

MUL2200 

MUL2210 

MUL2220 

MUL2230 

MUL2240 

HUL2250 

MUL2260 

HUL2270 

HUL2280 

MUL2290 

HUL2300 

MUL2310 

MUL2320 

MUL2330 

MUL2340 

MUL2350 

MUL 2 360 

MUL2370 

MUL23aO 

MUL2390 

MUL2400 

MUL2410 

MUL2420 

MUL2430 

MUL2440 

MUL2450 

MUL2460 

MUL2470 

MUL2480 

MUL2490 

MUL 25  00 

MI)L25lO 
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20  NSAW=NSAW*1 

CALL   SAW   IS&HVtPOINTtTOTAL*ISAH( IKAX) 
GO  TO  23 

21  NVENft=NVENR* 1 

CALL   VFNR    I VE NRV , P0| NT , TOT AL . I S AH . IH AX ) 

OQ  rn  23 

22  NSTUO=NSTUD* 1 

CALL    STUD    ISTUOV, OQINT, TOTAL. ISAK, IHAX ) 

23  CONTINUE 

CALl    PULP  (PULPV,P0INT,TI)T&L.ISAM,IMAX1 

LAST=2 
2^  CONTINUE 
25  CONTINUE 


C  BEGIN  ALL   PRODUCT  ROUTINES 


IFIJONREF.NE.3)    Gfl  TO  ^99 
26  J0N=2 

RFWIND  I 
ICARD=0 

IF  (IPUC.EQ.n  WRITE  (3,^31  PROHOT 
LAST=1 


-READ  AND  EVALUATE   DATA  FOR  SELECTED  PRODUCTS 


on   38   KX=Kf IN,KHAX,2 

LIH=KnUNT(KX) 

IF    (LIM.CO.O)   GO  TO  3B 

no  37   J=l  ,LIH 

no  27  1=1,10 

27  iNntn=o 

RE4D( I t K, I STRT, IDVYtOBH, IPLHT, ISwP, IPCR.  ICRK, IPFK, IFRK, IFS, ILS,  I  (  t 
IKNTT  (L.M.N  I  ,N=l  ,  31  .M=  l,'.|,L  =  l,3l,II2APlU,l  =  l.'i),I  TRhT  ,  L  E  AN  »  I  ROT 
IFIITRHT.LF.O)    WRirEI3.*.7U  ITBHT.IDf'.Y 
IF ( 1 TRHT.LE.O 1    GO   TO  37 

CALL    VOLUM  (KCUV8iKeFV9.KeFVl6,eAFV» 

 IF    RANDOM  NUMBER   GENERATOR  I RANF I  I )  )    IS  NOT  USED  TO  DETERMINE  SUB- 

 SAMPLE   FOR  ROT  R E DUC T 1  ON . 1 N SE RT   A  GO   TO  30  BRANCH  CARD  HERE.  R=0 

GO  TO  30 

IF    (AANFU).LT-PROftQT)  28.29 
?e  IR0T=2 
r.O  TO  30 

29  IRaT=0 

30  IF    IIBUG.EO.O)   GO  TO  32 
ICARO=ICARD*l 

WRITCll.'t'))  ICARD,ISTRT,PniNT(ISTRT),EF,K,OBH,IPLHT,ITRHT,ISWP, 
I  LEAN, ICRK, I  PGR » IFRK,  IPFK,  IFS.  ILS.  IROT , ( I 2^P (IH ) , I M  =  l , 4  I 
DO   3  1    L=l  .3 
on   31  M=l,4 

31  «K[TE    (3,51)    L.M, 1KN0T(L,M,N} ,N=1,3I 

32  IF    I (KX.eO.KMAX) . AND. ( J.EO.L IH) J  LAST=3 
IF    IKPriLE.EO.O)   GO  TO  33 

CALL    POLE    (PLHT2, POINT, TOTAL,  ISAM,  IMAX.PLHTl) 

33  IF    1K.S4W.EQ. 01    GO  TO  3* 

:-J  S  A  H  =  0 

CALL    SAW    I SAHV. POINT, TOTAL,  ISAM,  I  MAX) 
3^   IF    (RVENR.EQ.Ol   GO  TO  35 
NVFNR=0 

CALL   VENR    ( VE NRV , PO 1  NT , T OT AL ,  I  SAM,  I  MAX  1 

35  IF    (KSTUO.EO.Ol    GO  TO  36 
NSTUO=0 

CALL    STUD    ISTUnV, POINT, TOTAL,  ISAM,  IMAX » 

36  IF(KPULP.E0.01   GO  TO  60 

CALL    PULP    (PULPV, POINT, TOTAL,  ISAH,  IMAX  ) 
60    IFIKUUT.EQ.O)    GO   TH  62 

IFlKP.EO.l)    I PLHT=I PLHT*5 

WRl  TE  K,ISTRT  ,  I  DM  Y  ,  DBM  ,  I  PLHT  ,  IS  HP,  IPCR,  ICRK,  IPFK,  J  FRK,  IFS,  ILS, 

I  (  1  (KNOT  (L,M,M  ,N=l,3»  ,M=l,^  t  ,L  =  1,  3  1,  (  I  2^P  (  I  )  ,  I  =  1 ,  4  >  ,  I  TR  HT  ,  LE  A  N  , 
21ftOT,KP, ICFV, TBFV.EF, ( ( BF SW 1 1 , L 1 , L  =  I, 5  > , 1  =  1 , 4 ) , ( (BFVNI l.Ll,L=l,3), 
31  =  1        ,  (  IPFSTI  I  ,L)  ,L=1 .3)  ,  1  =  1        ,  ICFPL(  11.1  =  1.2) 

NREC=NHEC  ♦  I 
62  LAST=2 

37  CONTINUE 

38  CONTINUE 
ENDFILE  4 
HEHlNO  4 

IFIKOliT.EQ.O)   GO  TO  999 
WRlTh(3,53)  NREC 

WRI TE (3,54  1    K, JSTRT, lOMY, DBH, IPLHT, ISWP, IPCR. I CRK , I PFK , I FRK , I F S . 
IILS.M  (KN0T(L,M,NI,N=1,31  ,H=l,41,L  =  l,  31,  (  I24P(I  l,I  =  l,4),|TRHT, 
2LEAN, IROT,KP,TCFV,TflFv,EF , (  I RF SW (  1 , L I , L =  1  , 5  1  ,  I  =  1  ,  4  )  ,  I  ( PFVN(  I  ,L)  ,L  = 
3l,3),I=l,41,l(RFST(I,Ll,L=l.3),  1  =  1,4)  ,  ICFPH  11.1  =  1,21 

GO   TO  999 

39  CALL  EXIT 


^.0  FORMAT  (3A101 

41  FORMAT  11H0»SAMPLE  S I Z E* I 6 . ZX *L AHGER  THAN  ALLOWABLE,  SEE  DATA  STAT 
1EMENT*1 

42  FORMAT  (2I5,F5.3,12I5,F5.01 

43  FnRHAT(9I5,2X,Fq.41 

44  FORMAT  (16F5.0) 

45  FnRMAT(2XI2. 17, F5. 2, 1x12, 1X711, U12II,  2X1211,  2X1211, 2X4 II, 1X12, 
11X21  1 ) 

46  FORMAT  ( 1 H0*0 I  A MET E R •F6 . 2 , 2 X*EX C E EOS  MAXIMUM  ALLOWABLE,  DATA  FROM 
ITHIS    TREE    >^0T  USEn«) 

47  FORMAT  (IH  ^SAMPLE  NUHBE R • 1 6 , 2X«EXCEeDS  MAXIMUM  ALLOWABLE.  DATA  FR 
liT^   THESE   SAMPLES   NOT  USED*) 

^,71  FHR^'ATllH  *lRBt  HE  I  GH  T  •  I  4  ,  2  X  «E  XCE  EDS  MINIMUM  OR  MAXIMUM  FOR  VOL  CA 
ILCULATION,    DATA   FOR    TREE    in»18,«      NOT    USED*  1 

'.R  F^KM^T  I  IHl,  130(  4=*1  /  1H0«TREF  STATISTICS  PIECE  BY  P I  ECE  ^F  20.  5 , 3  X  •O 
IF    BUTT    LUGS  DISCARDED   PASEO  ON   RCiT    S  AMPL  F  «/ IHO ,  1  30  (  *  =•  1  1 

i.)  FOffMaTllH  ,130(*  =  «(/1H  'TREE  NUHBER*I6,«  SAMPLE  NUM8ER«I5,*  SUBP 
UlTS  IN  THIS  CLUSTCR»F5.0,«  VOL/ACRE  EXPANSION  F  ACTOR'F  9. 4 /•  TREE 
2    SIZE  DIAMETER    CLASS*I3,»      01 AME TFR*F6 . 2 , «      POL E   HE  I GHT* I  3 , *  T 

3-(EE  HEIGHT*|3/»  OFFECTS  SWEEP*)2,»  LFAN»I2,»  CR00K*I2,*  CROO 
4K  L0CATI0N*12,«  F0RK*I2,*  FORK  L OC A T 1  ON • I  2 , •  FIRE  SCAR*I2,»  LI 
5GHTNING  SCARCI2,"  R0T*I2/*  CLEAR  FACES  IN  FOURTH  P1ECF«4I2,/*  KNO 
6TS   PIECE  FACC  T,OTAL  LIVE  0EA0«1 

50  FORMAT    1 IH   ,  I  10,4191  " 

51  FORMAT    I  IH    ,  no  ,4  I-i  1 

53  FORMATdHO*  ENO-OF-Rl)N  TOTAL  RECORDS  OUTPUT  ON  TAPE  4  UNDER  OPTI 
im  2   =tl5,e  LAST   RECORD  PROCESSEIl  WAS  •! 

54  FORMAT (IHO, 13, 14, 18, F6.2,  13, 1X712,  1X1212,  1X1212,  1X1212,  1X412,  13, 
l?I2/lH0,XX!2,2F8.2,F1.4/lH0,2X,?0F6.Z/lHn,2x,20F6.2/lH0,2x,6F6.2) 

END 

SUBROUTINE    POLE    ( PLHT2 , PO I  NT , TOT AL , ISAM, !MAX,PLHTl  ) 

C  

r  ASSIGNS   ACCEPTABLE  POLES  TO  HEIGHT   AND  DIAMETER   CLASSES  AND 

C  ACCUMULATES  POLE  DATA 

C  

DIMENSION  PLHT21 ISAM. ISl .    POINT(ISAM),    TO T AL (  I S AM , ^  I  ,  PLHT1(15,20) 

C 

COMMON  KNOT (3,4,31,I?4P(41,!Pllll,INni  1 0 1 , KOUN T ( 50 ) , Y I  151  ,S(  15) , 
lPROB(fll ,PCUVfl( \ I  1 ,Pprvn(  lOI . PftFVl6( 5) , BFSW( 4, 5 l .ftFVN(4, 3)  , 
2tiFST(',,3),CFPL(2),NPK(m,KMIN,KHAX,NSAM,K,ISTRT,0RH,IPLHT,ISWP, 
3IPCR.ICHK,IPFK,IFHK,1FS,ILS,ITHHT,LEAN,  I  ROT , N S A W , N VE NR , N S TUD , JON , 


MUL2520 
MUL2530 
MUL2540 
MUL2550 
HUL2560 
HUL25  70 
HUL2580 
MUL2590 
MUL2600 
MUL2610 
MUL2620 
HUL2630 
HUL2640 
MUL2650 
MUL2660 
HUL2670 
HUL2680 
MUL2690 
HUL2700 
MUL2710 
MUL2720 
MUL2730 
HUL2740 
HUL2750 
MUL2760 
HUL2770 
MUL2780 
MUL2790 
MUL2800 
HUL2310 
MUL2820 
MUL2830 
MUL2840 
MUL  2850 
MUL2B60 
MUL2870 
MUL2e80 
HUL2890 
HUL2900 
HUL2910 
MUL2920 
HUL2930 
HUL2940 
MUL2950 
MUL2960 
MUL2970 
MUL2980 
MUL2990 
MUL 3000 
MUL3010 
MUL3020 
MUL3030 
MUL3040 
HUL3050 
MUL3060 
HUL3070 
HUL3080 
HUL3090 
HUL3100 
HUL3110 
MUL3120 
MUL3130 
HUL3140 
HUL3150 
MUL3160 
MUL3170 
HUL31B0 
MUL3190 
MUL3200 
MUL32iO 
MUL3220 
HUL3230 
MUL3240 
HUL3250 
HUL3260 
HUL3270 
MUL 3280 
MUL3290 
MUL3300 
MUL3310 
MUL3320 
HUL3330 
MUL3340 
MUL3350 
MUL3360 
HUL3370 
HUL3380 
HUL3390 
MUL3400 
HUL3410 
MUL3420 
MUL3430 
MUL3440 
MUL3450 
HUL3460 
MUL3470 
MUL3480 
MUL  3490 
MUL3500 
MUL3510 
MUL3520 
MUL3530 
MUL3540 
MUL3550 
MUL3560 
MUL3570 
MUL3560 
MUL3590 
MUL3600 
MUL3610 
MUL3620 
MUL3630 
MUL3640 
MUL3650 
HUL3660 
HUL3670 

POL  10 

POL  20 

POL  30 

POL  40 

POL  50 

POL  60 

P(tL  70 

POL  60 

PPL  90 

POL  100 

POL  110 


4LAST.TCFV,TBFV, I  BUG, I  CARD, PR  I SM, KP , I OMY. KVOL , EF 


KP  =  0 

r.n   TO    11.4.4  1,  LAST 

1  ltMN=10 

IF    (KHIN.GT.IO)  KMN=KMIN 

KMX=20 

IF    IKMAX.lT,20)  KMX=KMAX 

no   3  11=1.15 

UO  2  KK=KMN,KMX,2 

2  PLHTl  ( 1  I  ,KK)=0 
on  3  JJ=l,NSAM 

3  PLHT2 ( JJ, 1  I  1 =0 

4  IF    ( IPLHT, LE.O)   GO  TO  20 

IF    ( K.GT.20.0R.K.LT. 101    GO  TO  20 

IF1NK=( i IPLHT/8.01+0.91 

IF ( IFINK.GT. ITRHT 1    IFINK    =  ITRhT 


-GRADING  SECTION 

ITOP=0 
ICUT=0 

GO  TO   15,71 ,  JON 
DO  6  L=l . IFINK 
I F ( I ND 1 L  1  ) 

CONTINUE 


20 


7    IFdSHP)  8,8,20 

fl    IF ( I LS.GT.O.OR. IROT.GT, 1 1    GO  TO  20 

IFILEAN.GE.2)   GO  TO  20 

IFIICRK-1)  10,10,9 
9    IF(It>CR.GT.IFINKl    GO   TO  10 

IF(IPCR.NE.l.OR.lPCR.NE.IFINK)   GO  TO  20 

IF ( IPCR.EO. 1 1  ICUT=5 

IFII^'CR.Fg.IFINKl    ITGP  =  5 
10    IF  I  IFRK 1  14,14,12 
12    IF ( IPFK.GT. IF INKl    GO  TO  14 

IF ( IPFK. NE.l .OR. IPFK. NE. IFINK 1    GO   TO  20 

IFl  IPFK. FQ.l. AND. 1  CUT. FQ.O)    ICUT  =  5 

IF ( 1 PFK.EQ. IF  INK. AND. ITOP.eO.O)  IT0P  =  5 
14    IFIIFS.GT.O.ANO.lCUT.eO.O)  ICUT=5 

LIM=AMIN0I3, IFINK) 
DO   17  L=l,LIM 
00    17  M=l,4 

IF    (KNOT(L,H,2)-41  16,20,20 

16  IF    (KN0T(L,M,3)-4)  17,20,20 

17  CONTINUE 

IPLHT=I PLHT-ICUT-ITOP 
IFiNK=( ( I PLHT/8.0)+0.91 
IFdFINK.GT.ITRHTl    IFINK    =  ITRHT 
ZZ=CVOLtl.O, PRISM, EFI 


C  ADD   INTO  REGISTER 


IF ( IBUG.EO. 1 . AND.LAST.NE.3I    WR  ITE(  3, 28  1    I  CARD, IPLHT.ZZ, 

I PLHT= I PLHT/5 

IF    (  IPLHT. LT. 2)   GO  TO  20 

KP  =  1 

PLHTl  (  IPLHT,  K)  =  PLHTH  IPLHT,  K)*ZZ/POINT(  ISTRT) 
PLHT2 ( ISTRT, IPLHT 1=PLHT2( ISTRT, I PLHT ) +ZZ 
GO  TO   (  18,201 ,  JON 

18  00  19  L=1.10 

19  IN0(L)=1 

20  IF    (LAST.NE.3)   GO  TO  27 
WRITE    (3,291    IPROB( I  1  .  1  =  1,8) 
on  TO   (21,221,  JON 


C  PRODUCT   PRIORITY  HEADING 


Q  SINGLF   PRODUCT  HCAOING 


22  WRITE  (3,31) 

23  WRITE  (3,321 
DO  25  KK=KVN,KMX,2 
DO  24    I  1=2,15 

24  PLHTl  III  ,KK)=PLHTl(II  ,KK)/NSAM 
WRITE    (3,331    KK,  (PLHTII  II.KKl,  11  = 

25  CONTINUE 

no  26  11=2,15 

26  CALL   STDER    1 NSflM, Y I  I  I  I , S (  I  I  ) , PLHT 2 (  1 , 1 1  1 . PO INT  1 
WRITE    (3,34)    I Yl  11  1 , I  1  =  2,15  1 

MRI TE    (3,351    (  S (  1 1  1  ,  I  I  =2 , 1 5  1 

27  RETURN 


,15) 


28  FORMAT    IIH  •TREE   NUMBER* 1 6 , 1 2X ,* ACCEPTED  AS  POLE  CLASS*I3,2X, 
I   *WITH  PER   ACRE   WEIGHT  OF »F  1 0 . 3 , 2X , »0 1 AME TER  CLASS»I31 

29  FORMAT    ( IHI.RAIO) 

30  FORMAT    ( IHO, I  30 ( »=• ) / IHO , *POL ES    BASED  ON   PROUUCT  PRIORITY  CRITERIA 
I*) 

31  FORMAT    (1H0,130(*=*)/IH0,»PQLES  CONSIDERED  AS   SINGLE  PRODUCT*) 

32  FORMAT    ( IHO, 1 30 ( ♦=*) / 1 H0*0 1 AMETER*37X 'POL ES  PER   ACRE  BY  HEIGHT  CLA 
ISS*/1H0*  CLASS*16X 

2«10*5X*15«5X*20*5X*2S«5X*30«5X*35*5X*40*5X*45*5X*50*5X*55»5X*60* 
3   5X«65*5X*70«5X*759/i HO, 130 (*=*)) 

33  FORMAT    (IH    I  5  ,  1 3x , I 4F 7 . 2/ IH  1 

34  FORMAT    ( 1H0,130(*=*)/IH0»HE1GHT  CLASS   SUM  *.14F7.2) 

35  FORMAT    (1HO*STO   ERR   MF ANe6 X , 1 4F 7 . 2 / I HO ,  1 30  I *  =  * )  1 
END 

SUBROUTINE    SAW    ( VC NT , PO I  NT . TOT AL , I S AH . IM A X 1 

GRADES  LOGS  FOR   SAWMILL  -MATERIAL  USING  A  FIVE  GRADE  SYSTEM 
ACCUMULATES   SAWLOG   DATA   FOR    16FT   LOG  AND  HALF  LOG 

, IMAXl , POINT ( ISAM) 
L( ISAM, 4 1 

COMMON  KN0T(3,4,31 . I24P(4 1 ,  1 P  (  1  U  ,  I  NO (  1 0 1  .  K OUNT  (  50 ) . Y ( 1 5 ) . S ( 1 5 ) . 
1  PROP  181  ,PCUV8 (lll,PeFV8(10),PeFVl6l5),BFSW(4,51,BFVN(4,3l, 
2BFST('.,3),CFPL(21  ,  NPROD  ,  KM  I  N  ,  K  M  Ax  ,  NS  AM  ,  K  ,  ISTRT,  DBH,  IPLHT,  ISWP, 
3 IPCR,  ICRK,  IPFK,  IFRK, IFS, ILS, ITRHT, LEAN,  I R OT , N SA W , N VENR . NS TUD , JON . 
4LAST,TCFV,TBFV, lOUG, I C AR D , PR  I SM , K P , IOHY,KVOL,EF 

GO   TO    11,3,3)  .LAST 

1  KMN=12 
IF(KMIN.GT.12)  KMN=KMIN 
DO  2  KK=KMN,KMAX,2 

00  2  J=l,4 
DO  2    I  =  UNSAH 

2  VCNT( I , J,KK)=0.0 

3  IF (K.LT.  12  )   GO  TO  200 
IF(N5AW-n4,4,5 

4  DO  7  1=1,4 
no  7  J=l,^ 

7  BFSW{ ! , J)=0.0 


POL  120 
POL  130 
POL  140 
POL  150 
POL  160 
POL  170 
POL  180 
POL  190 
POL  200 
POL  210 
POL  220 
POL  230 
POL  240 
POL  250 
POL  260 
POL  270 
POL  280 
POL  290 
POL  300 
POL  310 
POL  320 
POL  330 
POL  340 
POL  350 
PCL  360 
POL  3  70 
POL  380 
POL  390 
POL  400 
POL  410 
POL  420 
PCL  430 
POL  440 
POL  450 
POL  460 
POL  470 
POL  480 
POL  490 
POL  SCO 
POL  510 
POL  520 
POL  530 
POL  540 
POL  550 
POL  560 
POL  570 
POL  580 
POL  590 
POL  600 
POL  610 
POL  620 
POL  630 
POL  640 
POL  650 
POL  660 
POL  670 
POL  680 
POL  690 
POL  700 
POL  710 
POL  720 
POL  730 
POL  740 
POL  750 
POL  760 
POL  770 
POL  780 
POL  790 
POL  800 
POL  810 
POL  820 
POL  830 
POL  840 
POL  850 
POL  860 
POL  870 
POL  880 
POL  890 
POL  900 
POL  910 
POL  920 
POL  930 
POL  940 
POL  950 
POL  960 
POL  970 
POL  990 
POL  990 
POL  1000 
POLlOlO 
POL  1020 
POL  1030 
P0L1040 
POL1050 
POL  1060 
PCL1070 
POL1080 
POL  1090 
POLUOO 
POLIUO 
POL  1120 
POLll 30 
SAW  10 
SAW  20 
SAW  30 
SAW  40 
SAW  50 
SAW  60 
SAW  70 
SAW  80 
SAW  90 
SAW  100 
SAW  110 
SAW  120 
SAW  130 
SAW  140 
SAW  ISO 
SAW  160 
SAW  170 
SAW  180 
SAW  190 
SAW  200 
SAW  210 
SAW  220 
SAW  230 
SAW  240 
SAW  250 
SAW  260 
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C    TEST    FOR    TWO   CIR   MORE   HALF    LOGS    IN  TREE 
C 

5  IF ( I TRHT-2 1200,6,6 

6  LOC=L 
LPCl=l 
LPC2=2 

IF ( I TRHT.GT. 2)   GO  TO  8 

KHT  =  0 

GO  TO  10 
8  KHT=l 
10  GO  TO    112, 12)  ,JON 
12  IF(INOILPCll*IND(LPC2)-l)l'.,A7,^7 

C 

C  GRADE   FIRST  LOG   (2  SECTIONS) 
C 

m  NP  =  0 

DO  22  L=l,2 

N0=0 

N1  =  0 

00  20  M=l,'. 

IF[KUOTiL,M, 1 )-ll 16,18,20 
16   NO=NO  +  2 

GO  TO  20 
18  Nl=Nl  t  1 
20  C0^4TINUE 

NP  =  NP*-N0+N1 
22  CONTINUE 

IF (NPI26,26,24 

24  GO  TO   (26,26,26,26,28, 28,28, 28, 28, 2R, 2B. TO, 30,30,32,32 ) ,NP 

C 

C    IGR=4     DENOTES   LOG  GRADE  5 
C 

26  IGR=4 

IRG  =  5 

GO  TO  34 
28  IGR=3 

I«G  =  3 

GO  TO  34 
30  IGR=2 

1«G  =  2 

GO  TO  34 
32  IGR=1 

1RG  =  1 

34  GO   TO    (36,40,36) , JON 
36  IFJNSAH~IGR)47,38,47 
38  1ND(LPC1)=1 
IN01LPC2)=l 

40  IF(KVOL-l)  903,903,900 

C 

C   ADJUST  GROSS   LOG  VOL   FOR  VISUAL   DEFECTS   FOR   16FT  LOG 
C 

900  PCNT=0-0 

IF(ILS.EQ.l)    PCNT=PCnT   +  .25 

IF(ILS.eQ.2)    PCNT=PCNT  +  .50 

IF( ISWP.E0.2)    PCNT=PCNT   t  .20 

IF { IFS.EQ.2. AND.LPCl.EQ.  1 )    PCNT  =  PCNT  +  .13 

I F ( I CRK.GT.O. AND. IPCR.EO.LPCl .OR. IPCR.EQ.LPC2  )  PCNT=PCNT  +  .25 
IF (I FRK.GT.O. AND. IPFK .FO. LPCl .Oft . IPFK,EQ,LPC2)    PCNr=PCNT    +  .25 

C 

C   ADD   INTO  REGISTERS 
C 

41  IFIPCNT-l.)  42,44,44 

42  VLM= (PBFV16I LOG)-( PBFV 16 ( LOG  )  *PCNT  )  ) 

43  VLM=CVOL(VLM, PRISM, EF) 
GO   TO  46 

44  VLM=0.0 
IF(LPCl.E0-LPC2)   GO  TO  45 
iriD(LPCl)=0 
IND1LPC2)=0 

r.D   TD  46 

45  INO(LPCU=0 

46  VCNT  (  I  STRT,  IGR.K  l=VCNTnSTRT,  IGR  ,K  )    +  VLM 
IF(IBUG.EQ.1.AND.LAST.NE.3)    GO   TO  904 

GO  TO  47 

C 

C   NO  GROSS  LOG  VOL   ADJUSTMENT   FOR  VISUAL   DEFECTS   16FT  LOG 
C 

903  Vt.M=PBFVl6(L0G) 

8FSW( IGR,LO& )=6FSH( IGR.LOG)   *  VLM 
GO  TO  43 

C 

904  GO   TO    (906,906) ,KVOL 

908    WRlTE(3,57n    I  C  AR  D  ,  L  OG ,  1  RG  ,  V  L  H ,  K  ,  PCNT 

571    FORMATdH    *TR6E    NU  M8E  R  *  I  6  ,  2  X  *      L  OG*  I  3  ,  2  X  «  ACC  E  P  T  ED   AS    SAHLOG  GRADE' 
1I3,2X*HITH   ADJUSTED   WT   VOL    OF *F 1 0 . 3 , 2 X*n  I  AMET E R   C L AS S » I  3 , 2 X *P R V  =' 
2F5.2) 
GO  TO  47 

906   WRITE (3,570)    I  CARD , LOG, I  KG, VLM, K 

570  FORMATdH  *THEE   NUM6E  R*  1 6  ,  2  X  «     LOG*  1 3 ,  2X  *  ACCEPT  ED  AS   SAWLOG  GRADE' 
1I3,2X*WITH  WEIGHTED  VOLUME  OF*F  10. 3, 2X*0 1 AMETER  CLASS*I3I 

C 

C   TEST  FOR   SECOND  LOG 
C 

47  CONTINUE 

GO  TO   148,48,98,98, 1 18, 1 18, 126, 126,200,200) ,LPC2 

48  1F(KHT)200,200,50 
50    IF( ITRHT-4152,82,82 

C 

C    1TRHT=3     NEXT   SECTION   IS  HALF  LOG  (8FT) 
C 

52  L0G=2 
LPC1=3 
LPC2=3 
KHT  =  0 

GO  TO    (54,56,54) , JON 
54    IFIIND(LPCl) )56,56,47 

C 

C   GRADE    A   HALF    LOG    SECTION  IRFT) 
C 

56  N0=0 
Nl  =  0 

no  62    M=l ,4 

IF (KNOT ( 3,M, 1 )-l )58,60,62 
58  NO=NO  +  2 

GO  TO  62 
60  Nl=Nl  +  I 
62  CONTINUE 

NP=NO+Nl 

NP=NP*2 

1  F (NPl 66  ,66 ,64 

64   GO  TO    (66,66,66,66,68,68,68,68,68,68,68, 70, 70, 70, 72, 72) ,NP 


DENOTES  GRADE  5 


66  IGR=4 
I  RG  =  5 
GO    TO  74 


SAW  270 
SAW  280 
SAW  290 
SAW  300 
SAW  310 
SAW  320 
SAW  330 
SAW  340 
SAW  350 
SAW  360 
SAW  370 
SAW  380 
SAW  390 
SAW  400 
SAW  410 
saw  420 
SAW  430 
SAW  440 
SAW  450 
SAW  460 
SAW  470 
SAW  480 
SAW  490 
SAW  500 
SAW  510 
SAW  520 
SAW  530 
SAW  540 
SAW  550 
SAW  560 
SAW  570 
SAW  580 
SAW  590 
SAW  600 
SAW  610 
SAW  620 
SAW  630 
SAW  640 
SAW  650 
SAW  660 
SAW  670 
SAW  680 
SAW  690 
SAW  700 
SAW  710 
SAW  720 
SAW  730 
SAW  740 
SAW  750 
SAW  760 
SAW  770 
SAW  780 
SAW  790 
SAW  800 
SAW  810 
SAW  820 
SAW  830 
SAW  840 
SAW  850 
SAW  860 
SAW  870 
SAW  880 
SAW  890 
SAW  900 
SAW  910 
SAW  920 
SAW  930 
SAW  940 
SAW  950 
SAW  960 
SAW  970 
SAW  980 
SAW  990 
SAWIOOO 
SAWlOlO 
SAW1020 
SAH1030 
S AW1040 
SAW1050 
SAW1060 
SAW1070 
SAW1080 
SAH1090 
SAWLIOO 
SAWl  1  10 
SAKl  120 
SAWl 130 
SAWl  140 
SAWU50 
SAWl 160 
SAWl 170 
SAW  1  1 80 
SAW1190 
SAW1200 
SAW1210 
SAW1220 
SAW1230 
SAW1240 
SAW1250 
SAH1260 
SAW1270 
SAW1280 
SAW1290 
SAW1300 
SAW1310 
SAW1320 
SAW1330 
SAW1340 
SAW1350 
SAW1360 
SAW1370 
SAW1380 
SAW  1 390 
SAW  1400 
SAH14 10 
SAW1420 
SAW1430 
SAW1440 
SAW1450 
SAW1460 
SAW1470 
SAW1480 
SAH1490 
SAW1500 
SAW1510 
S4W1520 
SAW1530 
SAW1540 


68  IGR=3 

IRG  =  3 

GO  TO  74 
7  0   IGR =2 

IRG  =  2 

GO  TO  74 
72  IGR=1 

IRG=1 

74   GO   TO    (76,80,76) , JON 
76  IF(NSAW-1GR)47,78,47 
78  1ND(LPC2)=1 
80  IF(KV0L-1)913,913,912 

C 

C   ADJUST  GROSS   LOG  VOL   FOR  VISUAL   DEFECTS   FOR   6FT  LOG 
C 

912  PCNr=0.0 

IF(ILS.EQ.l)  PCNT=PCNT  *  .25 
IF(ILS.EQ.21  PCNT=PCNT  +  .50 
I  F  1  I  SWP.E0.2  )    PCNT-PCNT   ■»-  .20 

[FtlCRK.GT.O. AND. 1  PGR . FQ.LPCl )  PCNT  =  PCNT  ♦  .50 
IF  1  IFRK.GT.O. AND, I PFK .EQ.LPCl )  PCNT  =  PCNT  ♦  .50 
GO   TO  41 

C 

C   NO  GROSS  LOG  VOL   ADJUSTMENT   FOR  VISUAL   DEFECTS  8FT  LOG 


C    nRHT  =  4      NEXT    SECTION    IS    FULL    LOG  (16FT) 

C 

82  L0G=2 
LPC1=3 
LPC2=4 

IF( I TRHT.GT. 4)    GO   TO  81 

KHT  =  0 
GO   TO  83 
81  KHT=l 

83  GO   TO    (84, 86, 84), JON 

84  I F(1N0(LPC1)+IND(LPC2)-1) 86,4  7,47 

C 

C   GRADE   SECOND   LOG  [16FT1 
C 

86  N0=0 
N1=0 
N2  =  0 

DO  92    M=l ,4 

IF(KN0T(3,M,1 )-l )83,90,92 

88  NO=NO  +  2 

GO  TO  92 
90  N1=N1  t  1 
92  CONTINUE 

DO  96  M=l,4 

IF(  I24P(M1  )94,94,96 
94   N2=N2    ♦  2 
96  CONTINUE 

NP=N0+N1+N2 

IF  INP) 26,26,24 

C 

C   TEST   FOR   SECTIONS  ABOVE   TWO  LOGS 
C 

98  IF(KHT)200,200,IOO 
100   IF ( ITRHT-6) 102, 108, 108 

102  L0&=3 
LPC1=5 
LPC2-5 

103  KHT=0 

GO   TO    1  104, 106, 104) , JON 

104  IF( !N0(LPC1) ) 106,106,47 
106  IGR=4 

IRG  =  5 
GO   TO  74 

108  L0G=3 
LPCl=5 
LPC2=6 

IF( ITRHT.GT.6)    GO  TO  110 

109  KHT=0 

GO  TU  112 

110  KHT=1 

112   GO   TO    I  1  14,  1  16, 1 14) , JON 

114  IF(IN01LPC1)+IN0(LPC2I-1I116,47,47 

116  IGR=4 

IR&  =  5 

GO   TO  34 

C 

118  IF(KHT)200,200,120 

120    IF ( I TRHT-8 ) 122, 124, 124 

122  L0G=4 

LPCU7 

LPC2=7 

GO  TO  103 
124  L0G=4 

LPC1=7 

LPC2=8 

IF( I TRHT.GT. 8)    GO  TO  110 
GO  TO  109 

C 

126  IF(KHT)200,200,128 

128  IF(ITRHT-10)130,132,132 

130  L0G=5 

LPCl=9 

LPC2-9 

GO   TO  103 
132  L0G=5 

LPC1=9 

LPC2=10 

GO  TO  109 

C 

200   1F(LAST.NE.3)   GO  TO  500 

WRI TE(3,580)    (PROBl  I  ) .  1  =  1 .8) 
GO  TO   1204,202) , JON 

C  

C    SINGLE    PRODUCT  HEADINGS 

C 

202   GO   TO  (914,916),KV0L 
914  WRITE(3,586) 

GO  TO  918 
916  WRITE{3,588) 
91R  WRITE(3,590) 

GO  TO  212 
204  lF(NSAW-4)206,208,206 
206  NSWW=NSAW 

GO  TO  210 
208  NSWW=5 

C 

C    PRODUCT   PRIORITY  HEADINGS 


210  GO  TO  (920,922) ,KVOL 
920   WRITEi3,596)  NSWW 


SAW1550 
SAW1560 
SAW1570 
SAW1580 
SAW1590 
SAW1600 
SAW1610 
SAW1620 
SAW1630 
SAW1640 
SAW1650 
5AW1660 
SAW1670 
SAW1680 
SAW1690 
SAW1700 
SAW1710 
SAW1720 
SAW1730 
SAW1740 
SAW1750 
SAW1760 
SAW1770 
SAW1780 
SAW1790 
SAW1800 
S AH1810 
SAW1820 
SAW1H30 
SAW1840 
SAW1850 
SAW1860 
SAW1870 
SAH1880 
SAW1890 
SAW  1900 
SAW1910 
SAW1920 
SAW1930 
SAW1940 
SAW1950 
SAW1960 
SAW1970 
SAW1980 
SAW1990 
SAW2000 
SAW2010 
SAW2020 
SAW2030 
SAH2040 
SAW2050 
SAW2060 
SAW2O70 
SAW2080 
SAW2090 
SAW2100 
SAW21 10 
SAW2120 
SAW2130 
S4H2140 
SAW2150 
SAW2160 
SAW2170 
SAW2180 
SAW2190 
SAW2200 
SAW2210 
SAW2220 
SAW2230 
SAW2240 
SAW2250 
SAW2260 
SAW2270 
SAW2280 
SAW2290 
SAW2300 
SAW2310 
SAW2320 
SAW2330 
SAW2340 
SAW2350 
SAW2360 
SAW2370 
SAW2380 
SAW2390 
SAW2400 
SAW2410 
SAW2420 
SAW2430 
SAW2440 
SAW2450 
SAW2460 
SAW2470 
SAW2480 
SAW2490 
SAW250O 
SAW2510 
SAW2520 
SAW2530 
SAW2540 
SAW2550 
SAW2560 
SAW25  70 
SAH2580 
SAW2590 
SAW2600 
SAW2610 
SAW2620 
SAH2630 
SAW2640 
SAW2650 
SAW2660 
SAW2670 
SAW2680 
SAW2690 
SAW2700 
SAW2710 
SAW2720 
SAW2730 
SAW2740 
SAW2750 
SAW2760 
SAW2770 
SAW2780 
SAW2790 
SAW2800 
SAW2810 
SAW2820 
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GO  TU  92*  SAH2a30 

922   WMITEn»59Bl    NSHW  SAH28*0 

92*.   HRITe(3,600)  SAW2850 

212  DO  21*    I=l,NSaM  SAW2B60 

on  21*.   J=l  ,*  SAH2870 

21<.    TOTAL  (  1  ,  J  )  =0.0  SAH2680 

00  226    KK=KMN,KMax,2  SAW2690 

no  216   J=l  ,1,  SAW2900 

21^   CALL    STOeB    I NSAH, Y IJ ) ,S1J I , VCNT ( 1 , J,KK 1 .POINT)  SAH2910 

GO  TO   (220, 21*1 1  .  JON  SAH2920 

218   HRlTei3,6lO)    KK , ( V ( [ ) , 1 = 1 , * )  Saw2930 

C  SAW29<.0 

C   SINGLF   PRODUCT  OUTPUT  SAH2950 

C  SAW2960 

WRITE!  1,6201    I  Sn  )  ,  1=        )  SAW2970 

GO  TO  222  SAH2980 

C  SAH2990 

C    PRODUCT   PBIQKtTY   OUTPUT  SAW3000 

C  SAW3010 

220   dRiTe(3,630l    KK.YINSAW)  S4H3020 

WRI TE ( 3,6*01    SINSAW)  SAH3030 

222  no  22*    J=l,*  SAH30*0 

10  22*    l  =  UNSAM  SAM3050 

22*    TOTALU  ,  J1  =  T0IAL  (  I  ,  J1*VCNT  (  I  ,  J.KK  »  SAH3060 

226  CONTINUE  SAH3070 

on   228    J=l.*  'iAW3080 

228   CALL    STOER    ( N5 AM , Y ( J  1 , S IJ 1 , TOT  At (  1  ,  J  1  , PO  INT  )  SAW3090 

GO   TO    (232  ,230)  , JON  SAH3100 

230  HRI  TE  (  3.6"iO)    I  Y  (  1  I  ,  1  =  1  ,    )  SAW3110 

WHl TE I  3 ,6601    1S(I1,I  =  1,*)  SAH3120 

GO  TO  500  Saw3l30 

232  MR|TEn,670)    YtNSAU)  S&U3140 

WR|TEI3,680»    SINSAH)  SAW3150 

500  RETURN  SAW3160 

C  Saw3170 

S80  FORMAT  (1H1,9A10(  SAW3180 
586  FOHMATl IHO,! 30) »=» 1/1H0*SAML0GS  CONSIDERED  AS  SINGLE   PRODUCT       GRO  SAW3190 

ISS   VOLUME    NOT    ADJUSTED  FOR   VISUAL    DEFFCTS»1  SAM3200 
588  FORMAT(IHO,130(*=*1/1HO*SAHLOGS  CONSIDERED  AS   SINGLE   PRODUCT       GRO  SAH3210 

ISS    VOLUME    ADJUSTFP   FOR    VISUAL    DEFECTS*!  SAH3220 
590  FORMATl  lHO,I301*  =  *l/lHO»[HAMETeRe32)f*ROARO  FOOT  VOLUMES   PER   ACRE    B  SAW3230 

lY  GRAOE«/IHO»    C  L  A  S  S*  2  5  X  e  1  *  1  7  X  •  2  «  I  7  X  *  3  *  1 7  X      «  /  I  HO  ,  1  30  I  •  =  •  ))  SAh32*0 
596  FORMAT  I  IHO,  1301  *=<- 1/1H0«GRADE*I3.  *    SAHLOr,   VOLUMES   BASED   ON   PRODUCT  SAW3250 
1    PRIORITY  CRITERIA      GROSS    VOLUME    NOT    ADJUSTED   FOR    VISUAL    DEFECTS*!  SAW3260 
598   FORMAT  I IHO , I  30 ( *=• 1 / lHn*GRAOE* 1 3, *    SAWLOG  VOLUMES    BASED  ON   PRODUCT  SAW3270 

1    PRIORITY   CRITERIA      GROSS   VOLUME    ADJUSTED   FOR    VISUAL   DEFECTS*)  SAW3280 
600  FORMAT! IHO, 1301 •=• 1/1H0«DIAMETER*17X*B0ARD  FOOT  VOLUMES   PER   ACRE    *  SAW3290 

1/IH    «    CLASS*/1M0.130I*=*1 )  SAW3300 

610   FORMAT    ( IHC, I*,6X, •MEAN*F20.2,3I F 1 6.2  1  1  SAW3310 

6?0  FORMAT    I  IM    ,10X*STD   E RR 'F  I  7 . 2 , 3  I F  1  8 . 2  1  )  Saw3320 

63PFORMai{lHO,I*,?*X«MEAN*F20.2  1  SAW3330 

6*0  FORMAT    (IH    ,28X«STD  £RR»FIT.2)  SAW3340 

650   FORMAT    I  I HO , 1 30  I •  =  « ) / 1 HO* T OT AL      M E AN *3 X , F 20 . 2 , 3  I F 1 8 , 2 )  )  SAW3350 

660  FORMAT    (lHO*STn  ERR    ME  AN* 2 X , F20 . 2 , 3 ( F 1 8 . 2  1  / I HO , I  30 ( *  =  • )  )  SAW3360 

(-70  FORMAT    UH0,130I*  =  *l/lHO*T0TAL*23X*MEaN*F20.2)  S  Art  3370 

680  FORMAT    tlH0,28X»ST0  ERR *F  1  7 . ? /  1  HO ,  I  30 1 •= »  )  )  SAW3380 

END  S4VJ3390 

SUBROUTINE  VENR    1 VCNT , POI NTt TOTAL . I SAH, 1  MAX )  VEN  10 

C   VEN  20 

C  GRADES  FOR  VENEER  POTENTIAL  USING  A  THREE  GRADE   SYSTEM  VEN  30 

C  ACCUMULATES  VENEER  VOLUMES  VEN  *0 

C   VEN  50 

DIMENSION  VCNT(ISAM,*,IMA)(),POINTI  I  S  AH  )  ,  TOTAL  (  I  SAM  ,  *  1  VEN  60 

C  VEN  70 
COMMON   KNOT I3,*,3).!2*P(*I,IP(III,IN0(  10  1  ,KOUNT( 50)  , Y (  1 5 ) , S t 1 5 1 ,        VEN  80 

lPRnB(eI.PCUV8(ll),PBFV8ri0),PBFV16l5),BFSH(*,51,BFVN(*t3),  VEN  90 

2BFST l*,3 1  .CFPL  12  I  ,  NPROD , K M 1 N , KMAX , NS AM, K ,  ISTRT.DBH, IPLHT,  ISWP.  VEN  100 
3IPCR,ICRK,lPFK,IFHK,IFStlLS,  lTRHr,LEAN,  I RO T , N SAH , N VENR , NS TUD , JON ,      VEN    1 10 

*LAST,TCFV,T6FV, I PUG , I C AR 0 , PR  I SM , K P , lOHY.KVOL.EF  VEN  120 

C   VFN  130 

GO  TO   (1.3,3),   LAST  VEN  1*0 

1  KMN=10  VEN  150 
IF  IKHIN.tiT.  10)  KMN=KMIN  VEN  160 
DO  2    KK=KMN,KMAX,2  VEN  170 

00  2  J=l,*  VEN  180 
DO  2    1=1 ,NSAM  VEN  190 

2  VCNT 1 1 , J , KK ) =0.0  VEN  200 

3  IF  IK.LE.IO)  GO  TO  21  VEN  210 
IF    (K-1*)    *,5.6  VEN  220 

*   LIM=1  VEN  230 

GO  TO  7  VEN  2*0 

5  LIM=2  VFN  250 
GO  TO  7  VEN  260 

6  LIM=3  VEN  270 

7  L|M=AM|N01LIM,ITRHT)  VEN  280 
|F(NVENR-l)8,a, 10  VEN  290 

8  DO  9  1  =  1  ,*  VEN  300 
DO  9   J=l»3  VEN  310 

9  OFVNI I , J)=0.0  VEN  320 
10  on  20  L=1,LIM  VEN  330 

GO   TO   (50, 52). JON  VEN  3*0 

50   IF ( INDIL ) )52,52,20  VEN  350 

C   VEN  360 

C  GRADING  SECTION  VEN  370 

C   VEN  380 

52    IGK=l  VEN  390 

DO  5*  M=l,*  VEN  *00 

IF(KNOTIL,M,3).GT.0.OR.KNOT(L,M,21.GT.21    I GR= AM AXO ( I GR . 2 )  VEN  *10 

IFlKfJ0T(L.M,3).GT.2.0R.KN0r(L,M,2).GT.2)    IGR=AMAX0(IGR,3)  VEN  *20 

IF(KNQT(L,M,3).GT.3)    K.Rs*  VEN  *30 

IFIKNOTIL.M,3).&T.3)    GO   TO  55  VEN  **0 

5*  CONTINUE  VEN  *50 

iF(K.Lr.l6)    IGR=AMAX0t IGR.2)  VEN  *60 

55  CONT INUE  VEN  *70 
GO  TO   (56, 60). JON  VEN  480 

56  IFINVENR-IGR)20,58,20  VEN  490 
58  IM0(L>=1  VEN  500 
60   IF (KVOL-l ) 102, 102, 100  VEN  510 

C   VEN  520 

C  ADJUST  GROSS   LOG  VQJ.  FOR  VISUAL  DEFECTS  VFN  530 

C   VEN  5*0 

100  PCNT=0.0  VEN  550 

IFIIROT.GT.l.ANO.L.EO.n    PCNT  =  PCNT  +   1.0  VEN  560 

1  F (  |C»K ) 6* ,6* ,62  VEN  570 
62    IF  1  ICRK.GE.2 . AND.L.EO. IPCR)    PCNT=PCNT    +    1.0  VEN  580 

IF (  ICRK.EO. l.ANO.L.EO. IPCR)    PCNT  =  PCNT   ♦   .5  VEN  590 

6'.    IF  (  I  FRK.GT.O.  AND.L.EO.  I  PFK)    PCNT  =  PCNT    ♦    .5  VCN  600 

IF  (  I  LS)68,69 ,66  VEN  610 

66    IF(ILS.CF.21    PCNT=PCNT    *    .5  VEN  620 

IF(ILS,EQ.ll    PCNI=PCNT    ♦    .25  VEN  630 

68    IF(IFS.GE.2.AN0.L.FQ.l)    PCNr=PCNT    *    .25  VEN  6*0 

C   VEN  650 

t  ADD   INTO  REGISTERS  VEN  660 

C   VEN  670 

IF(PCNT-1.)70,72,72  VEN  680 

7U    VLM= 1PBFV8 (L l-( PBFVai L 1*PCNT ) I  VFN  690 

71    VLM=CVOL ( VLM.PR ISM.EF 1  VEN  700 

GO    TQ    7*  VFN  710 


72  VLM=0.0  VEN  720 

1ND(LI=0  VEN  730 

7*   VCNT I  I STRT ,IGR,K)  =  VCNT ( ISTRT,  IGR .K )    ♦    VLM  VFN  7*0 

If ( I BUG.EO. 1 . AND.LAST.NE. 3)    GO  TU    10*  VEN  750 

GO   TO  20  VEN  760 

C  VEN  770 

C   NO  GROSS  LOG  VO'    ADJUSTMENT   FOR  VISUAL   DEFECTS  VEN  780 

C  VEN  790 

C    102    IFdROT.GT.l.AND.L.EC.n    GO  TO  99  vEN  800 

C             IF (  ICRK.GT. I .AND.L.EO.  IPCR)    GO   TO  99  VKN  810 

C             1 F (  IFRK.GT.O. AND.L.EO.  IPFK)    GO   TO   99  VEN  820 

C            IFdLS.GT.O)    GO   TO  99  vEN  830 

C            IFI IFS.GE, 2. AND.L.EO. I)    GO  TO  99  vEN  8*0 

102   VLH=PPFVB(L1  VEN  850 

fiFVN( IGR,Ll=eFVN( IGR.L )    ♦  VLM  vEN  860 

C'O    TO    71  VEN  870 

C  VEN  880 

C    SET   IND   SELECTOR  OFF     PIECE  FAILED  VISUAL    DEFECTS   TEST  VEN  890 

C  VEN  900 

C            IN0(L)=0  VEN  910 

C           GO  TO  20  VEN  920 

C  VEN  930 

10*  GO   TO    ( 108,1061 ,KVOL  VEN  9*0 

106  WRIT£(3,36l)    t C ARO , L , I GR , VL M , K , PCNT  VEN  950 

GO   TO   20  VEN  960 

108   HRITE13,36)    I C ARD , L , I GR , VLM , K  VEN  970 

20  CONTINUE  VEN  980 

21  IF  (LAST.NF.3)  GO  TO  35  VEN  990 
HRI TE  (3.37)  ( PROB ( I  1 , I  =  1 , 8  I  VENIOOO 
GO  TO   (23,221,   JON  VENlOlO 

C   VEN1020 

C  SINSLE   PRODUCT  HEADING  VtNl030 

C   VEN10*0 

22  GO  TO  lll0,ll2)tKV0L  VEN1050 
no   WRITE13,376)  VEN1060 

GO  TO   11*  VEN1070 

112   WRITE(3,3781  VENIOBO 

U*    HRITE(3,381  VEN1090 

GO  TO  2*  VENllOO 

C   VENlllO 

C  PRODUCT    PRIORITY  HEADING  VFN1120 

C   V6N1130 

23  GO  TO  (116,118).KV0L  VEN11*0 
116  WRITE(3,386)    NVENR  VFN1150 

GO  TO  120  VEN1160 

118   WRITei3,388)    NVENR  VEN1170 

120  HR!T£(3,39)  VEN1180 

2*  DO  25   J=l,*  VEN1190 

00  25    I=l,NSAM  VEN1200 

25  T0TAL( I , Jl =0.0  VEN1210 
DO  31  KK  =  KMN,K,MAX,2  VEN1220 
DO  26   J=l,*  VEN1230 

26  CALL  STOER  1 NS AM, Y ( J  I  , S ( J  1  , VCNT 1  1 , J, KK ) , PO INT )  VEN12*0 
GO  TO   (28,27),   JON  VEN1250 

C   VEN1260 

C  SINGLE   PRODUCT  OUTPUT  VEN1270 

C   VEN1280 

27  WRITE  (3,*01  KK , I Y 1  I  1 , 1  =  1 , * )  VEN1290 
WRITE  I3,*ll  (S(Il,I=l,*l  VEN1300 
GO  TO  29  VEN1310 

C   VEN1320 

C  PRODUCT  PRIORITY  OUTPUT  VEN1330 

C   VEN13*0 

28  WRITE  13,*21  KK.YINVENRl  VEN1350 
WRITE    (3,*31    S(NVENR)  VeN1360 

29  DO  30  J=l,*  VEN1370 
on   30    1=1 ,NSAM  VEN1380 

30  TOTAL! I , J1=T0TAL( 1 , J)+VCNT(  I  ,  J,KK  )  VEN1390 

31  CONTINUE  VEN1*00 
00  32   J  =  l  ,*  VEN1*10 

32  CALL  STDER  1  NSAM ,  Y  I  J  I  ,  S  1  J  )  ,  TOT  AH  1 ,  J  )  ,  PO  I  NT  1  VEN1420 
GO  TO    (3*, 33),    JON  VEN1*30 

33  WRITE  (3,**)  (Yd), 1  =  1,*)  VEN1**0 
WRITE  (3,*5)  (S(I),I=1,*)  VEN1*50 
GO   TO  35  VEN1*60 

3*  WRITE    (3,*6)    Y(NVENR)  VEN1*70 

WRITE    (3,*7)    S(NVENR)  VENl*80 

35  RETURN  VEN1*90 
C   V6N1500 

36  FORMAT  ( IH  »TREE  NUMBER* 1 6 , 2X *P I  EC E * I  3 , 2 X*ACCEP T ED  AS  VENEER  GRADE  VEN1510 
1»I3,2X*W1TH  WEIGHTED   VOLUME   OF*F 1 0 . 3 , 2X • D I  AM E TER   CLASS*I3)  VEN1520 

361    FORMATdH   •TREE    NUMB  ER  *  I  6  ,  2  X  ep  I  EC  F  *  I  3  ,  2X  *  ACC  FP  T  EO   AS   VENEER    GRADE*  VEN1530 

II3,2X*WITH  ADJUSTED  HT  VOL  OF*F 1 0 . 3, 2X*D ! AMET ER  CL ASS* I  3, 2X»PRV  =«  VEN15*0 

2F5.2)  VEN1550 

37  FORMAT  dHl,8A10)  VEN1560 
376  FORMAT! IHO, 130(*=*1/1H0*VENFFR  CONSIDERED  AS   A  SINGLE  PRODUCT     GflO  VENISTO 

ISS  VOLUME   NOT   ADJUSTED  FOR  VISUAL   DEFECTS*)  VEN1580 
378  FORMATdH0,13O(*=*l/lH0*VENEER  CONSIDERED  AS  A   SINGLE   PRODUCT     GRO  VEN1590 

ISS  VOLUME   ADJUSTED  FOR  VISUAL   DEFECTS*)  VEN1600 

38  FORMAT( 1H0,130I*=«)/ lHn*0IAHFTER*23X«6QAR0  FOOT  VOLUMES  PER  ACRE  B  VEN1610 
lY  GRA0E*/lH0»   CL A S S *? 5 X * 1  *  1 7X « 2 •  1 7X* 3*  1  7X *** / IHO ,  1  30 ( •  =  • )  )  VEN1620 

386  FnRMAT( IHO, 1301 *=« 1/1H0*GRA0E*! 3,*  VENEER   BASED  ON  PRODUCT  PRIORIT  VEN1630 

lY  CRITERIA      GROSS    VOLUME    NOT    ADJUSTED    FOR   VISUAL    DEFECTS*!  VEN16*0 

388  FORMAT( 1H0,130(*=*I/1H0*GRA0E*I3,*  VENEER   BASED  ON  PRODUCT  PRIORIT  VEN1650 

lY  CRITERIA     GROSS  VOLUME   ADJUSTED  FOR   VISUAL   DEFECTS*)  VEN1660 

39  FOrtMATdHO,l30(*=*l/lHO*OIAMETER*17X*80ARO  FOOT  VOLUME  PER  ACRE  •  VEN1670 
1/lH   *   CLASS*/IH0,130(*=*1 )  VEN1680 

*0  FORMAT    (  IHO, I*  .6X*MFAN*F20.2,3(F18.2I 1  VEN1690 

*l    FORMAT    IIH    ,10X*ST0  ERR * F 1 7 . 2 , 3  I F  1  8 . 2  I  I  VEN1700 

*2    FORMAT    ( IHO, I*,2*X«HEAN*F20.2)  VEN1710 

*3   FORMAT    (IH    ,28X*STD  ERR*F17.2)  VEN1720 

**    FORMAT    ( IHO , 1 30  I *=  » 1 /IHO , *TOTAL      MEAN*3X, F20.2,3IF18.2))  VEN1730 

*5  FORMAT    (IH0*STD  ERR   M E AN* 2 X , F 20 . 2 , 3 ( F 1 8 . 2  1  / IHO , 1 30 ( •  =  • 1 1  VEN17*0 

*6  FORMAT    (lH0,13n(*=*l/ 1 HO »T0 F AL • 2 3X *ME AN* F 20 . 2 )  VEN  1750 

*7   FORMAT    (IH0,2BX*STD  E RR *F 1 7 . 2 / 1 HO , 1 30  I •  =  • ) )  VEN1760 

END  VEN1770 

SUBROUTINE   STUO   (VCNT , POI NT . TOT AL , ISAM, IMAX )  STU  10 

C   STU  20 

C  GRADES  FOR  STUO  POTENTIAL  USING  A  THREE  GRADE   SYSTEM  STU  30 

C  ACCUMULATES  STUD  VOLUMES  STU  *0 

C   STU  50 

DIMENSION  VCNTI ISAM,*. 9) .POINT (  ISAM)  , TOTAL(  ISAM,*!  STU  60 

C  STU  70 
COMMON  KNOT( 3,*,3 1 , 1 2 *P ( * ) ,  I P (  1 1)  , I  NO (  I  0 1 ,KOUNT (501.Y(151,S(I5!,        STU  80 

IPROB 181 ,PCUV8 (in .PBFV8 (10),P6FV16(51,BFSW(*,S1,BFVN(*,31.  STU  90 

2BFST(*,31,CFPL(2) , NPRQD , K M I N , KM A X , NS AM , K , I S TR T , OB H , I P L HT , ISWP,  STU  100 
31  PC R,ICRK,IPFK,IFRK,IFS.ILS,ITKHT, LEAN, IROT,N SAW, NVENR, NSTUD, JON,      STU  110 

*LAST,TCFV, TBFV,ieUG,ICARD,PRISM,KP, IDMY,KVOL,EF  STU  120 

C   STU  130 

GO  TO   11,3,3).   LAST  STU  1*0 

1  KMN=6  STU  150 
IF  )KMIN.GT.6l  KMN=KMlN  STU  160 
KMX=22  STU  170 
IF  (KMAX.LT.221  KMX=KHAX  STU  180 
00  2  KK=1,9  STU  190 
DO  2  J= I ,*  STU  200 
DO  2    I=1,NSAM  STU  210 

2  VCNT ( I , J,KK)=0.0  STU  220 
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KHN=(KHN-'.)/2  STU  230 

KMX=(KMX-^)/2  STU  Z^tO 

3  IF  (K.LT.6)  GO  TO  22  STU  250 
KK=(K-4)/2  STU  260 
IF    (K-22)    5,4,22  STU  270 

4  LIM=3  STU  280 
tin   TO   12  STU  290 

5  IF    {K-20)    7,6,6  STU  300 

6  LIH=2  STU  310 
GO  TO  12  STU  320 

7  LIM=1  STU  330 
12   UHT=AM1N0( ITRHT,3)  STU  3^0 

IF    (LIK.GT-LHT)   GO  TO  22  STU  350 

IF(NSTUO-lia,8,10  STU  360 

fi   DO  9    1=1,4  STU  370 

DO   9   J=l  ,3  STU  380 

9  8FST 1 1 , J)=0.0  STU  390 

10  DO  21    L=LIM,LHT  STU  400 

GO   TO    150,52)tJ0N  STU  410 

50   IF(IND(L1)52,52,21  STU  420 

52  Nl=0  STU  430 

C   STU  440 

C  GRADING    SECTION                                                                                    .  STU  450 

C   STU  460 

IGK^l  STU  470 

DO   54    M=l,4  STU  480 

N1=NH-KN0T  (  L  ,M,  I  I  STU  490 

IF (KNOT (L,M,3 ) . GT . 1 . OR . KNOT ( L , M , 2 ) .GT.2 1    1GR=AMAX0( IGR,2)  STU  500 

IF  (KNOT  (L,  M,3  1  .GT.2.0R.KN0T(  L.M,  2  )  .GT.2)    IGR=AMAJ<0(  I  GR  ,  3  1  STU  510 

IF (KN0TIL,M,3 ) .GT.2.0R.KN0T( L,M,2 ) .GT.3)    IGR=4  STU  520 

1FIKN0T(L,M,3).GT.2.0R.KN0T(L,M,21.GT.3)    GO  TO  55  STU  530 

54  CONTINUE  STU  540 
IFIN1.GT.16)  IGR=AMAXO( IGR,2)  STU  550 
IF(N1.GT.32I    IGR=AMAXO( IGR,3)  STU  560 

55  CONTINUE  STU  570 
GO   TO    156, 60), JON  STU  580 

56  IF(NSTU0-IGR)21  ,58,21  <.  STU  590 
58  INDH.)  =  1  STU  600 
60    IF (KVOL-1  )  102,  102,100  STU  610 

C   STU  620 

C  ADJUST  GROSS   LOG  VOL   FOR  VISUAL   DEFECTS  STU  630 

C   STU  640 

100   PCNT=0.0  STU  650 

IF(IROT.GT.l.ANO.L.EQ.l)    PCNT=PCNT   +    1.0  STU  660 

IF(ICRK)64,64,62  STU  670 

62    IF ( ICRK.GE.2 . AND.L.EO. IPCR)   PCNT=PCNT   +   1.0  STU  680 

IF(  ICRK.EQ.l  .AND.L.EO.  IPCR)    PCNT  =  PCNT   ■•■   .5  STU  690 

64    IF  1 IFRK.GT.O. AND.L.EO. I PFK)    PCNT=PCNT    +    .5  STU  700 

IF ( ILS )6n,6B,66  STU  710 

66    IF(1LS.GE.2)    PCNT=PCNT   +    .5  STU  720 

IFIlLS-eO.l)    PCNT=PCNT    +    .25  STU  730 

68    IFl 1FS.GE.2. AND.L.FQ. I )    PCNT=PCNT    ♦    .25  STU  740 

C   STU  750 

C  ADD    INTO  REGISTERS  STU  760 

C   STU  770 

IF  (PCNT-1. )70,72,72  STU  780 

70  VLM=  1PBFV8 (L )-( PRFV8 (L )*PCNT )  )  STU  790 

71  VLM^CVOLIVLM, PRISM, EFl  STU  800 
GO   TO   74  STU  810 

72  VLM=0.0  STU  820 
IND(L)=0  STU  830 

74   VCNT( I STRT  ,  IGR,KK)=VCNT(  ISTRT , ICR.KK I    +   VLM  STU  840 

IF ( I6UG.EQ.I. AND.LAST.NE. 3)    GO  TO    104  STU  850 

GO  TO  21  STU  860 

C  STU  870 

C   NO  GROSS  LOG  VOL   ADJUSTMENT   FOR  VISUAL  DEFECTS  STU  880 

C  STU  890 

C    102   IF( IROT.GT. 1. AND.L.EO. 1)    GO  TO  99  STU  900 

C            IFl ICRK.GT.l .AND. L.EQ. IPCR)    GO  TO   99  STU  910 

C             I F I IFRK.GT.O. AND. L.ED. IPFK)    GO  TO   99  STU  920 

C             IF(ILS.GT.O)    GO   TO  99  STU  930 

C             IFI IFS.GE. 2. AND.L.EO. I)    GO   TO   99  STU  940 

102   VLM=PeFV81L)  STU  950 

BFST( IGR.L  )  =  BFST  I  IGR,L)    *   VLM  STU  960 

GO  TO  71  STU  970 

C  STU  980 

C   SET   IND  SELECTOR   OFF     PIECE  FAILED  VISUAL   DEFECTS  TEST  STU  990 

C  STUIOOO 

C      99    IND(L)=0  STUlOlO 

C           GO  TO  21  STU1020 

C  STU1030 

104  GO    TO    ( 108,106) ,KVOL  STU1040 

106   WR1TE13,371)    I C AR 0 , L , 1 GR , VL M , K , PCNT  STU1050 

GO  TO  21  STU1060 

108   WRITE(3,37)    1 C AR D , L , I GR , V L M , K  STU1070 

21  CONTINUE  STU1080 

22  IF  (LAST. NE. 31  GO  TO  36  STU1090 
HRI TE  13,38)  ( PROB( I  1,1  =  1,8)  STUllOO 
GO  TO   (24  ,23),   JON  STUIUO 

C   STU1120 

C  SINGLE   PRODUCT  HEADING  STU1130 

C   STU1140 

23  GO  TO  (  110, 112) ,KVOL  STU1150 
110  WRI  TE  (3,3861  STUH60 

GO    TO   114  STU1170 

112   HRIT£(3,389I  STU1180 

114  WR1TE(3,39)  STU1190 

GO  TO  25  STU1200 

C   STU1210 

C  PRODUCT   PRIORITY  HEADING  STU1220 

C   STU1230 

24  GO  TO  1 116,118) ,KVOL  STU1240 
lift  WRITE(3,396)    NSTUD  STU1250 

GO  TO  120  STU1260 

118  WRITE(3,398)    NSTUD  STU1270 

120  WRITE(3,40)  STU1280 

25  no  26  J=l,4  STU1290 
DO  26    1=1,NSAM  STU1300 

26  TOTAL( I , J)=0.0  STU1310 
DO  32  KK=KHN,KMX  STU1320 
KX  =  KK<'2*4  STU1330 
00  27    J=l,4  STU1340 

27  CALL  STDER  ( NSAM, Y ( J ) , S 1  J  ) , VCNT (  1 , J , KK ) , PO ! NT )  STU1350 
CO  TO   (29,26),    JON  STU1360 

C   STU1370 

C  SINGLE   PRODUCT  OUTPUT  STU1380 

C   STU1390 

20    WRITE    (3,41)    KX  ,  ( Y (  1  1  ,  I ^  1  ,  4  )  STU1400 

WRITE    (3,42)    (Sill, 1  =  1, 4)  5TU1410 

GO   TO  30  STU1420 

C   STU1430 

C  PRODUCT    PRIORITY  OUTPUT  STU1440 

C   STU1450 

29  WRITE  (3,43)  KX,Y(NSTUD)  STU1460 
WRITE    13,44)    SINSTUD)  STU1470 

30  DO  31  J  =  l  ,4  STU1480 
00  31    1=1, NSAM  STU1490 

31  TOTAL ( I  , J  1 = TOTAL ( I , J) *VCNT(  I  ,  J.KK  )  STU 1500 


32  CONTINUE  STU1510 
DO   33   J=l,4  STU1520 

33  CALL  STDER  ( NSAH  ,  Y ( J ) , S ( J ) , TOT AL ( 1 , J ) , PO I  NT )  STU1530 
GO  TO    (35,341,    JON  STU1540 

34  WRITE  (3,45)  (Yd), 1  =  1, 4)  STU1550 
WRITE  (3,46)  (S(l),I=l,4)  STU1560 
GO  TO  36  STU1570 

35  WRITE  (3,47)  Y(NSTUD)  STU1580 
WRITE    (3,48)    S(NSTUD)  STU1590 

36  RETURN  STU1600 
C   STU1610 

37  FORMAT  (IH  *TREE  NUMBER* 1 6, 2X*P I ECE*I 3 , 2X*ACCEPTED  AS  STUD  GRADE  STU1620 
1*I3,2X«W1TH  WEIGHTED  VOLUME  OF*F 10 .3 , 2X*0 1 AMETER  CLASS*!3)  STU1630 

371  FORMATdH  *TREE  NUMBE  R  *  1 6 ,  2X«P  I  EC  E«  1 3 ,  2X  «  ACCE  PTED  AS  STUD  GRADE*  STU1640 
1I3,2X*WITH   ADJUSTED  WT   VOL    OF *F 1 0 . 3 , 2 X *D I  AM E T ER   CL AS S » 1 3 , 2X«PR V  =*  STU1650 

2F5.2)  STU1660 

3B   FORMAT    tlHl,8Al0)  STU1670 

386  F0RMAT(1H0,130(*=»1/1H0*STU0S  CONSIDERED  AS   SINGLE  PRODUCT     GROSS  STU1680 

IVOLUME    NOT   ADJUSTED   FOP   VISUAL    DEFECTS*)  STU1690 

388   FnRMAT(lH0,130(*-»l/lH0*STUDS  CONSIDERED  AS   SINGLE   PRODUCT     GROSS  STU1700 

LVOLUME    ADJUSTED   FOR   VISUAL    OEFECTSO  STU1710 

39  FORMAT ( IHO , 1 30( *=* 1 / 1H0*0IAMETER*23X*60ARD  FOOT  VOLUMES  PER  ACRE  B  STU1720 
lY  GRADE*/ IHO*  CLASS*25x*l*l7X*2*17Xt3*l7x*-i,*/lH0,  130 (»  =  •})  STU1730 

396   FORMAT(  1H0,130(*  =  *)/IH0*GRADE*I3,*   STUDS    BASED  ON   PRODUCT   PRIORITY  STU1740 

1    CRITERIA     GROSS    VOLUME    NOT    ADJUSTED  FOR    VISUAL    DEFECTS*)  STU1750 

39R   FORMAT  I  IHO, 130( *  =  • )/IH0*GRADE*I3, *   STUDS    BASED  ON   PRODUCT    PRIORITY  STU1760 

I    CRITERIA      GROSS    VOLUME    ADJUSTED   FOR   VISUAL    DEFECTS*)  STU1770 

40  FORMAT ( IHO, 130( *=*1/IH0*D1AMETER*17X*B0ARD  FOOT  VOLUME  PER  ACRE  *  STU1780 
I/IH   *    CLASS*/ IHO, 130( *  =  *)  )  STU1790 

41  FORMAT    ( IHO, I4,6X*HEANtF20.2,3( F18.2  )  )  STU1800 

42  FORMAT    (IH    ,  10X*ST0  E RP *F 1 7 . 2 , 3 ( F 1 8 . 2 )  )  STU1810 

43  FORMAT    (  1 HO , I  4 , 2 4 X *ME AN*F 20 . 2 )  STU1820 

44  FORMAT  I IH  ,28X*STD  ERR*F17.2)  STU1830 
4  5   FORMAT    (IHO, 1 30  I  *  =  * ) / 1 HO *TOT AL *2 X *ME AN * , 3 X , F 20 . 2 , 3 ( F I  8 , 2 )  )  STU  1840 

46  FORMAT    (1HO*STO   ERR    ME AN * 2 X , F 20 . 2 , 3  I F 1 B . 2 ) / 1 HO , I  30 ( *  =  * ) )  STU1850 

47  FORMAT    ( 1  HO , 1 30 ( *= * ) / I HO *T0 T AL «2 3 X *ME AN* F 20 . 2 )  STU1860 

48  FORMAT  (IH0,28X*ST0  ERR*F 1 7 . 2/ LHO ,  1 30 1 *  =  *  1  )  STU1870 
EUQ  STU1880 

SUBROUTINE    PULP    ( PULP V , PO I NT , TOT AL , I S AM , I M AX )  PUL  10 

C   PUL  20 

C  ACCUMULATES  CUBIC  FOOT  VOLUME  FOR   PIECES  NOT  USED  AS   OTHER  PRODUCT  PUL  30 

C   PUL  40 

DIMENSION   PULPVI 1  SAM) , POINT ( ISAM) , TOT AL( ISAM, 4)  PUL  50 

C  PUL  60 
COMMON  K NOT ( 3,4, 3) , I24P(4) ,  I P 1  1 1 ) ,  I  NO (  10  )  , K OUN T ( 50  )  , Y (  151  , S(  15) ,        PUL  70 

1PR0B(81  ,PCUV8( in  ,P8FV«(  101  ,PBFV16( 5 ) , DFS H ( 4 , 5 ) , BF VN ( 4 , 3 1 ,  PUL  80 

2BFST(4,3I,CFPL(2)  , NPROD, KM  I N, KMAX, NS AM, K ,  I S TR T , 06 H , I P L HT , 1 SWP ,  PUL  90 
31PCR, ICRK, IPFK, IFRK, IFS, ILS, ITRHT, LEAN, IROT ,N SAW , N VENR , NS TUD , JON,      PUL  100 

4LAST,TCFV,TBFV,IBUG,ICARD,PRISM,KP,IDMY,KV0L,EF  PUL  110 

C   PUL  120 

GO  TO   (I, 20, 20), LAST  PUL  130 

1    DO    15    1=1, NSAM  PUL  140 

15  PULPVl I )=0.0  PUL  150 

20   NT=0  PUL  160 

KTOP=ITRHT   ♦    1  PUL  170 

00  25   1=1,2  PUL  180 
25  CFPH  I  )=0.0  PUL  190 

DO   30    1  =  1  .ITRHT  PUL  200 

30   NT  =  NT   +    INDl  I  1  PUL  210 

IF( NT)  35,35,85  PUL  220 

C   PUL  230 

C  TREE   EXCLUDED  FROM  PRIOR   PRODUCT  CONSIDERATION  PUL  240 

C  MAKE   GROSS  VOLUME  ADJUSTMENTS   FDR   VISUAL   DEFECTS  PUL  250 

C   PUL  260 

35    IF(KV0L-l)60,60,37  PUL  270 

37  DO   50    1=1, ITRHT  PUL  280 

PCNT=0.0  PUL  290 

IF(1LS)42,42,40  PUL  300 

40   PCNT=PCNT   ♦    1.0  PUL  310 

42    IF ( IFS.GT.O. AND. I .EO.  1)    PCNT  =  PCNT   +   .5  PUL  320 

IF(PCNT-L. 144, 46, 46  PUL    3  30 

44   VLM= ( PCUVR( I  )-( PCUV8(  I  ) *PCNT  )  )  PUL  340 

GO   TO  48  PUL  350 

46  VLM=0.0  PUL  360 

48    VLM=CVOL(VLM, PRISM, EFl  PUL  370 

PULPV(ISTRT)=PULPV(ISTRT1    *  VLM  PUL  380 

IF( IBUG.EQ. I. AND.LAST.NE. 3)    WRITE(3,10l)    1  CARD. I , VLH, K , PCNT  PUL  390 

50  CONTINUE  PUL  400 

t   PUL  410 

C  00  TOP  PIECE   FOR  ADJUSTED  GROSS   VOLUME  PUL  420 

C   PUL  430 

51  IF  (  ILS)52,52.54  PUL  440 

52  VLM=PCUV8 ( KTOP)  PUL  450 
PCNT=0.0  PUL  460 
GO  TO   56  PUL  470 

54   VLM=0.0  PUL  480 

PCNT=1.0  PUL  490 

56  VLM=CVOL(VLM, PRISM, EFl  PUL  500 

PULPV ( I STR T ) =PULPV ( I STRT 1    +   VLH  PUL  510 

I F( IBUG. EO. 1 .AND.LAST.NE. 31    WRITE(3,1031    I C AR D , VL M , K , PCN T  PUL  520 

58    IF(LAST.NE.3)    GO  TO  900  PUL  530 

PUL  540 

C  PUL  550 

C  00  TREE   FOR  GROSS  VOLUMES— NO  ADJUSTMENT   FOR  VISUAL  DEFECTS  PUL  560 

C  PUL  570 

60  00  62    1=1, ITRHT  PUL  580 

VLM=PCUV8 ( I )  PUL  590 

CFPL(21=CFPL(2)    ♦    VLM  PUL  600 

VLM=CV0L(VLM, PRISM, EF)  PUL  610 

PULPV( I STRTl =PULPV( ISTRTl    *   VLM  PUL  620 

IF ( IBUG.EQ. 1. AND.LAST.NE. 31    WRITE(3,10)    I C ARO , I , VL M , K  PUL  630 

62    CONTINUE  PUL  640 

C  PUL  650 

C  on   TOP    PIECE    FOR  GROSS    VOLUME — NO   ADJUSTMENT   FOR    VISUAL    DEFECTS  PUL  660 

C  PUL  670 

64   VLM=PCUV8(KT0P)  PUL  680 

CFPL ( 1 )=CFPL( 1 )    +   VLM  PUL  690 

VLM=CVOLlVLM, PRISM, EF)  PUL  700 

PULPV( ISTRT)=PULPV( ISTRT)    ♦   VLM  PUL  710 

IF( IBUG.EQ. 1. AND.LAST.NE. 3)    WRITE(3,91    ICARDrVLH,K  PUL  720 

GO   TO  58  PUL  730 

C   PUL  740 

C  STEM   SECTIONS  EXCLUDED  FROM  PRIOR  PRODUCT  CONSIDERATION  PLUS  TOP       PUL  750 

C  MAKE  GROSS  VOLUME   ADJUSTMENTS  FOR  VISUAL   DEFECTS  PUL  760 

C   PUL  770 

85    IF(KVOL-l )98,98,87  PUL  780 

87   00  96    1  =  1, ITRHT  PUL  790 

1  F(  IND(  I  )  .EO. I )  Gn  TO  96  PUL  800 
1  F  (  I LS) 90,90,92  PUL  810 

90  VLM=PCUV8(I1  PUL  820 

PCNT=0.0  PUL  630 

GO  TO  94  PUL  840 

92   VLM=0,0  PUL  850 

PCNT=1.0  PUL  860 

94    VLM=CVOL(VLH, PRISM, EFl  PUL  870 

PULPV( ISTRT) =PULPV( I STRT )    *   VLM  PUL  880 

IFIlBUG.EQ.l. AND.LAST.NE. 3)    WRITf(3,101)    I C ARD ,  I , VL M , K , PC N T  PUL  890 

96   CONTINUE  PUL  900 
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-DO  STEM  FOR  GROSS   VOLUMES — NO  ADJUSTMENT   FOR  VISUAL  DEFECTS 


18  on    100    1=1,1 TRHT 

IF  (INOm  .EQ.  U   GO  TO  100 
VLM=PCUV8( I ) 
CFPl.(21=CFPH2l    ♦  VLM 
VLM=CV0L(VLM,PR1SM,EF) 
PULPVI ISTRT)=PULPV(1STRT)  + 
IFIIBUG.EO.l. AND.L AST . NE. 3 ) 

100  CONTINUE 
GO    TO  6^ 

no   WRITEO.II)    I  PROBt  U  ,  1  =  1  ,8) 
GO   TO    (6,SI ,  JON 


-ALL   PRODUCTS  HFADING 


5  Gfl  TO  (102,  10^*)  ,KVOL 
102   WRl IE  <3  . 12) 

GO  ro  7 

10*.  HR!TE13,12H 
GO  TO  7 

6  GO  TO  ri06,  100) .KVQL 
106   WRI Tt 13 ,1 3) 

CO   TO  7 
108   WRI  TEO.  131  ) 

7  on  a  I=l,NSAH 

S    TDTALt I , I (=0.0 

CALL    STDER  (NSAH.Yd), 
WftI TE    (3,1^)    VI  1 1  ,S(  1 ) 
900  RETURN 


S  (  1  1  .PIILPVI  U, POINT  ) 


<>  FORMATIIH    *TREe   NUMBE  R«  I  6  ,  2  X«  T  OP   P  1  ECE*  .  I X*  AC  C  EP  T  ED  AS  PULP*13X»HI 
ITH  WEIGHTED   VOLUME    OF »F 1 0 . 3 , 2x *0 1 AME T ER  CLASS*I3) 

10  FORMAT    (  IH  "TREE    NUMBE R *  I  6 . 2X * P I E C E « I  3 , 2 X 4 AC C E P TE D  AS  PULP*13X*WIT 
IH   WEIGHTED   VOLUME   OF ^F I  0 . 3 , 2 X * D  I  AM E T FR  CLASS«I3) 

101  FORMATdH  'TREE  NUMBE  R  »  I  6  ,  2  X  "  P  I  EC  E*  I  3  ,  2X  *  ACC  E  PTEO  AS  PULP*13X»WITH 
1   ADJUSTED  WT   VOL   OF »F  1 0. 3 , 2 X *D I  AM E T ER  CL ASS* I  3 , 2X *PR V  =*F5.2) 

103  FORHATllH  STREE  NUMBER • I  6 , 2 X  * T OP  P I  EC E * , 1 X * AC CE P T ED  AS  PULP*13X»HI 
ITH  ADJUSTED  WT  VOL   OF *F 1 0 . 3 , 2 X *D I AME T ER  CL AS S« 1 3 , 2 X *PR V  =*F5.2) 

11  FORMAT  (IH1,8A10) 

12  FORMAT! IHO, 1301 •  =  «  )/lH0*PULP  CONSIDERED  AS  SINGLE  PRODUCT  GROSS  V 
lOLUME    NOT    ADJUSTED   FOR  VISUAL    OE EEC T S » / I  HO,  1 30 ( *  =  *  )  ) 

121  FORMAT! 1HO,130(»=«) /IH0»PULP  CONSIDERED  AS  SINGLE  PRODUCT  GROSS  V 
lOLUME    ADJUSTED   FOR   VISUAL    DE F EC T S ♦ /  1  HO ,  1  30 ! *  =  *  )) 

13  FnkMAT(IH0.l3O!»=*)/lHn«PULP   CONSIDERED   AS  RESIDUAL   PRODUCT  ONLY 
IGROSS   VOLUME    NOT    ADJUSTED   FOR   VISUAL    DE F EC T S» / 1 HO , I  3D  I «  =  *  )  ) 

131    FORMAT  I IHO, 130) *  =  »)/lH0»PULP   CONSIDERED  AS  RESIDUAL    PRODUCT  ONLY 
IGROSS  VOLUME   ADJUSTED  FOR   VISUAL   DE F EC T S « / I HO , 1 30 ! •=« ) ) 
M  FORMAT    IIH  *TOTAL   PULP  VOLUME    IN    CUBIC    FE E T »F 1 8 . 2/ 1 HO* S T AND AR 0  ERR 
iaR*F35.2/lHD,130I*  =  «l  ) 
END 

SUBROUTI NE    VOLUM    ( KCUV8 , KBFVB  »KBFV  l6. BAFV ) 

THIS  ROUTINE   CALCULATES  GROSS   CUFT   AND  BRET  TREE  VOL  USING  HGT   IN  16FT 
LOGS   AND  DBH  AND  ALSO  CALCULATES   A  VOL. /ACRE   EXPANSION  FACTOR  FOR  EACH 
SAMPLE   TREE   USING  06H  AND  THE    INPUT   BAF  VALUE  OF  THE  ANGLE  GAUGE  USED. 
THE   CALCULATED  GROSS  VOLUMES   ARE   ALLOCATED  AMONG  STEM  SECTIONS  OF  THE 
THEE    BY   SELECTED   PERCENTAGE  TABLES. 


DIMENSION   KCUVR(65) , K BF V8 I  55 ) , K BE V  16 ( 29 ) 

COMMON  KNOT  !  3, -V,  3)  .  IZ^Pf*),  IP(  11  )  ,  INDI  1 0  )  .  KOUNT  I  50  )  ,Y!15J.S!15), 
IPROBIB)  ,PCUVR  I  I  I  1  tPBFV8(  101 ,PBFV16( 5) , BF S W t A , 5 ) . B F VN ! ^ , 3 )  , 
2BFSTU,3).CFPL12I  ,  NPROD,  KM  I  N,  KHAX  »  NSAM,  K  ,  I  S  TR  T  ,  OB  H,  I  P  L  HT  ,  ISWP, 
31 PCR,  ICRK,  IPFK, IFRK, IFS, I L S , I TRHT , L E A N ,  I R OT , NS AW , N VE NR , NS TUO , JON , 
4LAST,TCFV,TBFV, 1  BUG , I C ARD , PR  I SH , KP ,  IDMY,KVOL. EF 

A  =  0.0 
EF=0.0 
no  5   N=l , 11 
5  PCUV8IN)=0.0 
DO   10  N=l,lO 
10  P6FV8(N)=0.0 
00    15  N=l,5 

15  PBFV16IN)=0.D 

ASSUME   BRKNG    POINT   OF    19.0    INCH   DBH  FOR    BLACKJACK   AND  OLD  GROWTH  TREES 
TOT  TREE   VOL  EQUATIONS  FROM  MY E R S , C . A , , US F S  RES.   PAPER  RM-2»1963. 
CONVERT   HALF-LOG  HTIITRHT)    TO   LOG  HT(16FT  SECTIONS) 

HT16=FL0AT(ITRHT(/2. 
0SQR=DBH»*2 
IFIPRISM)  13,13,14 

13  01 V=OSOft*0. 005454 
EF=BAFV/DIV 

GO   TO  16 

14  EF=PR!SM 

16  A=DSQR*HTI6 
IFlOBH-19.0)  17,18,18 


PUL  910 
PUL  920 
PUL  930 
PUL  940 
PUL  950 
PUL  960 
PUL  970 
PUL  980 
PUL  990 
PULIOOO 
PULlOlO 
PUL  1020 
PUL1030 
PUL1040 
PUH050 
PUL  1060 
PUL1070 
PUL1080 
PUL 1090 
PULl 100 
PUL 11 10 
PUL  1120 
PUL 1130 
PUL  1140 
PUL 1 150 
PULl 160 
PUL1170 
PUL  1 130 
PUL1190 
PUL12D0 
PUL1210 
PUL1220 
PUL1230 
PUL1240 
PUH250 
PUL1260 
PUL1270 
PUL1280 
PUL1290 
PUL 1300 
PUL1310 
PUL1320 
PUL  1330 
PUL1340 
PUL1350 
PUL 1360 
PUL137D 
PUL1380 
PUL1390 
PUL 14D0 
PUL1410 
PUL1420 
PUL1430 
PUL144D 
PUL 1450 
PUL1460 
VOL  10 
VOL  20 
VOL  30 
VOL  40 
VOL  50 
VOL  60 
VOL  70 
VOL  80 
VOL  90 
VOL  100 


17  IF1A.GT.800. )    TC FV=0 . 04 42 04 • A   +  8.266 
IF!A.CT.800. )    TBFV=0.300081*A  -  52.090112 
IF  !  A.LE.800.  )    TCFV  =  0.046*A   +  6.8 

IF  I  A.LE.800. 1    T6FV  =  0.224793*A   +  8.1656 
GO   TO  19 

18  IF  lA.GT. 1000. )    TCFV  =  0.O45736eA   +  10.857212 
1FIA.GT.1130. )    TBFV=0.326427*A  -  62.962331 
IF! A.LE, 1000. )    TCFV=0.050666*A  ♦  5.8668 
IFIA.LE.  1130.  )    T8FV=0.2757e4*A  -  5.09125 


C  CUBIC  FT  VOL  8Y  HALF-LOG  THIS   TREE  +  TOP 
C 

19  KT0P  =  1TRHT  *■  1 
IF( ITRHT-1 120,20,25 

20  ISKP=1 
GO  TO  35 

25  ISKP=0 

DO   30  N=2. ITRHT 
30   ISKP=ISKP   +  N 

rSKP=ISKP  ♦  1 
35  00  40  N=l,KTOP 

PCUV8!N)  =  (TCFV9KCUV8!  ISKP)  1/100. 

|SKP=ISKP  +  I 
40  CONTINUE 

C 

C    BOARD   FT   VOL    BY  HALF-LOG  THIS  TREE 
C 

KTOP=ITRHT  -  1 

IF!ITRhT-1145,45,50 
45  1SKP=1 

GO   TO  60 
50  1SKP=0 

DO  55  N=l  ,KTOP 
55    I SKP=ISKP  ♦  N 

ISKP=!SKP   *  I 
60  DO  65  N=l, ITRHT 

PBFV8(N)=ITBFV*KBFV8(  ISKP  1  1/100. 

ISKP=1SKP   +  1 
65  CONTINUE 

C 

C   BOARD  FT  VOL  BY  SAW  L0G(16FT)   THIS  TREE 
C 

LP=HT16  -  0.5 

IFIHT16-1.  1100,70,75 
70  LS=1 

GO  TO  90 
75  N=HT16 

X=N 

D=HT16  -  X 

IF!0)80.P0,85 
80  LS=N**2 

GO  TO  90 
85  LS=N»*2   ♦  N 
90  KTOP=LP  +  I 

00  95  N=l,KTOP 

PBFV16IN1=(TBFV*KBFV16(LS) 1/100. 
LS=LS  +  1 
95  CONTINUE 
100  RETURN 
END 

FUNCTION  CVOL   ( VBL , PR  I S H , EF )  CVO  10 

c   cvn  20 

C  ROUTINE   TO  CONVERT  VOLUME  TO  VOLUME/ACRE   FOR  PLOTLESS   CRUISING  CVO  30 

C  AND    TO  CONVERT   VOLUME    OF   FIXED  SIZED   PLOTS   ACCORDING   TO  A   CONSTANT  CVO  40 

C   CVO  5C 

IFIPRISM. NE. 0.0)    CV0L=V6L*PRISM  CVO  60 

IFIPRISM. NE. 0.0)    GO  TO  5  CVO  70 

CV0L=V8L«EF  CVO  BO 

5  RETURN  CVO  90 

END  CVO  100 


VOL  400 
VOL  410 
VOL  420 
VOL  430 
VOL  440 
VOL  450 
VOL  460 
VOL  470 
VOL  480 
VOL  490 
VOL  500 
VOL  510 
VOL  520 
VOL  530 
VOL  540 
VOL  550 
VOL  560 
VOL  5  70 
VOL  580 
VOL  590 
VOL  600 
VOL  610 
VOL  620 
VOL  630 
VOL  640 
VOL  650 
VOL  660 
VOL  670 
VOL  680 
VOL  690 
VOL  700 
VOL  710 
VOL  720 
VOL  730 
VOL  740 
VOL  750 
VOL  760 
VOL  770 
VOL  780 
VOL  790 
VOL  800 
VOL  810 
VOL  820 
VOL  830 
VOL  840 
VOL  850 
VOL  860 
VOL  870 
VOL  880 
VOL  890 
VOL  900 
VOL  910 
VOL  920 
VOL  930 
VOL  940 
VOL  950 
VOL  960 
VOL  970 
VOL  980 
VOL  990 


SUBROUTINE    STDER    ( NS AM, Y , S , AMTRX , POI NT  1 

C  CALCULATES    STANDARD  ERROR   OF  THE    MEAN    FOR    CLUSTER    SAMPLING  WITH 

C  UNEQUAL    NUMBER  OF    SUBPLOTS   PER  CLUSTER 

DIMENSION  AMTRXll),  POINTIl) 

C  NSAM=NUM6ER  OF    SAMPLE    CL U S T £ R S ! I NPUT ) 

C  Y  =  OUTPUT  MEAN 

C  S=OUTPUT   STANDARD  ERROR  OF  THE  MEAN 

C  AMTRX=INPUT  DATA  VECTOR  OF   SIZE  NSAM 

C  POINT=NUMBER   OF   SAMPLE  POINTS  FOR   NSAM  CLUSTERS 

P=FLOAT!NSAH) 

Y=0.0 

G=0.0 

DO    I  NC=1,NSAM 

IF    IPOINT(NC ) .EO.0.0)   GO  TO  I 
DIV=AMTRXINC )/POINT!NCl 
V=Y*DIV*DIV 
G=G+D1 V 
1  CONTINUE 

X6=IY-(G*G1/P1/IP-L.J 

S=S0RTIXB/P1 

Y=G/P 

RETURN 

END 


STD  10 
STD  20 
STD  30 
STD  40 
STD  50 
STD  60 
STD  70 
STD  80 
STD  90 
STD  100 
STD  110 
STD  120 
STD  130 
STD  140 
STD  150 
STD  160 
STD  170 
STD  180 
STD  190 
STD  200 
STD  210 
STD  220 
STD  230 
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KULTI PKOCUCT    IMVENTOKY   PKOGRAM      TEST   DATA      SAMPLH  OUTPUT 


POLES  COMSIDEREP  AS   SINGLE  PRODUCT 


DIAr'ETEK 
CLASS 


10 


15 


20 


POLES  PER  ACRE  BY  HEIGHT  CLASS 
25  30  35  40  45 


50 


55 


60 


65 


70 


75 


12 
lA 
16 
18 
20 


0.00  0.00 

0.00  0.00 

0.00  0.00 

0.00  0.00 

0.00  0.00 


0.00 
0.00 
0.00 
0.00 
0.00 


.5B 
.  52 

0.00 
0.00 
0.00 


.61 
.47 

0.00 
0.00 
0.00 


0.00 
.49 
0.00 
0.00 
0.00 


0.00 
0.00 
0.00 
0.00 
0.00 


0.00 
0.00 
0.00 
0.  00 
0.00 


0.00 
0.00 
0.00 
0.00 
0.00 


0.00 
0.00 
0.00 
0.00 
0.00 


0.00 
0.00 
0.00 
0.00 
0.00 


0.00 
0.00 
0.00 
0.00 
0.00 


0.00  0.00 

0.00  0.00 

0.00  0.00 

0.00  0.00 

0.00  0.00 


HEIGHT  CLASS  SUM 
STD   ERR  MEAM 


0.00  0.00  0.00  1.09  1.08 
0.00        0.00        0.00  .55  1.08 


.49  0.00  0.00  0.00  0.00  0.00  0.00  0.00  0.00 
.49       0.00        0.00        0.00       0.00        0.00       0.00       0.00  0.00 


MULTIPRODUCT   INVENTORY   PROGRAM     TEST   DATA     SAMPLE  OUTPUT 


SAWLOGS  CONSIDERED  AS   SINGLE   PRODUCT       GROSS  VOLUME  NOT  ADJUSTED  FOR  VISUAL  DEFECTS 


DIAMETER 
CLASS 


BOARD  FOOT  VOLUMES   PER  ACRE  BY  GRADE 


12 


MEAN 
STD  ERR 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


292.60 
40.  11 


MEAN 
STD  ERR 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


300.44 
73.05 


MEAN 
STD  ERR 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


53.38 
53.38 


MEAN 
STD  ERR 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
•  0.00 


20 


MEAN 
STD  ERR 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


60.08 
60.  08 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


80.42 
80.42 


102.95 
93.  11 


MEAN 
STD  ERR 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


MEAN 
STD  ERR 


0.00 
0.00 


0.00 
0.00 


49.58 
49.58 


63.  10 
63.10 


TOTAL  MEAN 
STD  ERR  MEAN 


0.00 
0.00 


0.00 
0.00 


109.66 
109.66 


892.89 
67.47 


MULTIPRODUCT   INVENTORY  PROGRAM     TEST  DATA     SAMPLE  OUTPUT 


PULP  CONSIDERED  AS  SINGLE   PRODUCT     GROSS  VOLUME  NOT  ADJUSTED  FOR  VISUAL  DEFECTS 


TOTAL  PULP  VOLUME  IN  CUBIC  FEET 
STANDARD  ERROR 


234.43 
13.63 
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MULTIPRODUCT    IMVEMTORY   PROGRAM     TEST   DATA     SAMPLE  OUTPUT 


VENEER  COMSIOERtD  AS   A  SPICLE   PRODUCT     GROSS  VOLUME  NOT  ADJUSTED  FO^  VISUAL  DEFECTS 


DIAMETER 
CLASS 


BOARD  FOOT  VOLUMES  PER  ACRE  BY  GRADE 
I  7  3 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MFAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


0.00 
0.00 


0.  00 
0.00 


0.00 
0.00 


0.00 
0.00 


0. 00 

0.00 


0.00 
0.00 


0.00 

0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 

0.00 


37.  07 


94.52 
31.46 


0.  00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 

0.00 


43.43 
43.43 


103.00 
54.40 


73.61 
18.20 


121 .03 
58.31 


42.70 
42.  70 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


14.48 
14.48 


0.00 
0.00 


0.00 
0.00 


0.  00 

o.on 


0.00 
0.  00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


0.00 

0.  00 


12.02 
12.02 


0.00 
0.00 


MEAN 
STD  ERR 


0.00 
0.00 


65.  36 
65.36 


TOTAL  MEAN 
STD  ERK  MEAN 


0.00 
0.00 


380.42 
165.55 


251.82 
113.56 


12.02 
12.02 


MULTIPROnUCT    INVENTORY   PROGRAM      TEST    DATA      SAMPLE  OUTPUT 


STUDS  CONSIDERED  AS  SINGLE  PRODUCT     GROSS  VOLUME  NOT   ADJUSTED  FOR  VISUAL  DEFECTS 


DIAMETER 
CLASS 


BOARD  FOOT  VOLUMES   PER  ACRE   BY  GRADE 


MEAN 
STD  ERR 


42.06 
42.06 


115.54 
63.70 


32.93 
32.93 


93.98 
8.76 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


MEAN 
STD  ERR 


.  78 
.78 


76.67 
45.17 


0.00 
0.00 


0.00 
0.00 


17.85 
17.85 


0.00 
0.00 


0.00 
0.00 


159.55 
75.95 


42.70 
42.70 


0.00 
0.  00 


MFAN 
STD  ERR 


0.00 
0.00 


0.00 
0.00 


0.00 
0.00 


MEAN 
STD  ERR 


0.00 
0.00 


0.00 
0.00 


TOTAL  MEAN 
STD  ERR  MLAN 


51. B3 
38.  12 


192.21 
108.25 


50.  78 
50.78 


296.23 
111.14 


—  CSU,  Ft.  Collins 
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